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SECTION  1 

INTRODUCTION  AND  BACKGROUND 

1.1  INTRODUCTION 

The  MILL  RACE  High  Explosive  Test  was  conducted  on  September  16, 

1981,  at  White  Sands  Missile  Range,  New  Mexico.  Sponsored  by  Field  Com¬ 
mand  of  the  Defense  Nuclear  Agency,  it  consisted  of  a  number  of  experiments 
and/or  engineering  tests  by  both  Government  Agencies  and  private  contractors. 
Teledyne  Brown  Engineering  (TBE)  was  prime  contractor  for  the  Contrast 
Attenuation  Technique  (CAT)  experiment  which  is  the  subject  of  this  Report. 

Previous  work  by  TBE  under  DNA  sponsorship  led  to  the  development 
of  the  CAT  to  deduce  dust  densities  from  photographic  films.  The  objective 
of  the  MILL  RACE  experiment  was  primarily  to  validate  the  CAT  methodology  which 
could  then  be  used  to  obtain  quantitative  estimates  of  dust  densities  produced 
during  nuclear  weapons  testing  at  the  Nevada  Test  Site  (NTS) .  It  was  not 
the  intent  of  the  MILL  RACE  CAT  test  to  obtain  and  use  the  data  for  scaling 
from  non-nuclear  to  nuclear  cases. 

Because  of  a  number  of  questions  concerning  the  absolute  accuracy 
of  the  CAT  in  deducing  dust  densities  from  NTS  films,  this  test  was  de¬ 
signed  to  reproduce  as  closely  as  possible  the  experimental  apparatus 
and  physical  conditions  found  for  many  of  the  NTS  tests  in  a  carefully  con¬ 
trolled  and  specified  parametric  manner.  Data  obtained  from  this  experi¬ 
ment  would  be  used  to  determine  both  the  accuracy  and  sensitivity  of  the 
CAT  and  thus,  hopefully,  provide  DNA  and  its  contractors  with  a  reliable 
method  for  using  the  NTS  film  data  to  support  current  phenomenology  studies. 

Figure  (1-1)  depicts  one  of  two  similar  test  stations  fielded  by  TBE 
for  the  MILL  RACE  Test.  The  two  camera  enclosures  contained  a  total  of 
three  cameras — two  in  the  enclosure  nearest  GZ  and  one  in  the  other  enclo¬ 
sure.  All  three  cameras  had  fields  of  view  containing  all  of  the  targets; the 
camera  nearest  GZ  (camera  A)  was  used  as  the  reference  line  of  sight  so  that  the 
six  targets  were  symmetrically  spaced  within  its  field  of  view.  The  six 
targets  were  designed  to  provide  equal  angular  dimensions  when  viewed  by 
camera  A  of  approximately  4°  x  4°,  with  adequate  spacing  between  targets 
to  measure  the  adjoining  background  radiance. 

Also  shown  in  Figure  1-1  are  several  dust  sampling  and  measuring 
instruments.  Although  not  part  of  the  TBE-fielded  equipment,  these 


Fig. (1-1)  Layout  of  the  CAT  validation  experiment  at  MILL  RACE.  Shown  are  two  camera  enclosures 
(labeled  1  and  2)  and  six  photo-targets.  Camera  enclosure  1  container  two  cameras — a  GSAP  nearest 
GZ  and  a  Flight  Research.  The  optical  centerline  passed  through  the  GSAP  camera  and  through  the 
center  of  the  photo-targets  located  at  17*  and  90’.  The  two  units  denoted  with  asterisks  on  the 
center  were  PMS  dust  sampling  and  measurement  instruments  used  for  real  time  dust  density  and 
particle  size  distribution  determination. 


measurements  were  conducted  by  SA1  as  an  integral  component  of  the  total 
validation  of  the  CAT.  In  addition,  at  the  station  nearest  GZ  (  -20  psi)  , 
SRII  fielded  a  microwave  transmission  experiment  within  the  field  of  view 
of  our  cameras  to  obtain  attenuation  vs  dust  density  data  during  MILL 
RACE;  this  also  was  part  of  the  CAT  validation  experiment. 

Two  test  stations  were  fielded  to  insure  that  sufficient 
dust  would  be  generated  at  one  or  both  sites  for  the  CAT  validation  test. 
Selection  of  the  two  locations  was  governed  by  several  considerations: 
first,  the  validation  experiment  was  predicated  on  obtaining  a  large 
enough  dust  cloud  to  encompass  the  cameras,  targets  and  all  of  the  sur¬ 
rounding  area  and  at  a  high  enough  density  to  cause  total  obscuration 
of  targets  as  close  as  ten  feet  from  the  cameras;  second,  the  experiment 
had  to  survive  the  effects  of  the  detonation  and  continue  to  operate 
thereafter.  After  studying  data  from  several  previous  H.  E.  tests,  it 
became  clear  that  a  site  at  20  psi  would  satisfy  the  first  requirement 
while  allowing  us  to  meet  the  second  requirement  with  reasonable  confi¬ 
dence.  Because  we  were  not  absolutely  convinced  that  survival  was 
guaranteed  at  20  psi,  we  decided  to  place  a  second  test  station  at  10  psi 
which  was  far  enough  away  from  GZ  to  guarantee  survival  but  not  clearly 
within  the  heavy  dust  cloud  as  predicted  from  previous  tesLs.  It  turned 
out  that  usable  data  was  obtained  from  both  stations  although  the  20  psi 
station  suffered  some  damage,  albeit  not  irreparable  as  far  as  validating 
the  CAT. 


1.2  BACKGROUND 

Nuclear  produced  dust  has  stymied  the  best  efforts  of  many 
researchers  to  accurately  predict  its  formation,  motion,  intensity  and 
dispersion  under  various  conditions  of  yield,  height-of-burst ,  surface 
conditions,  etc.  For  many  military  systems,  both  tactical  and  strategic, 
the  presence  of  large  dust  concentrations  can  severly  degrade  the  perfor¬ 
mance  of  radar,  IR  sensors,  missiles,  etc.  Recently  there  has  been  a 
fairly  substantial  renewed  effort  to  build  computer  models  that  accurately 
predict  the  dust  environment  produced  during  nuclear  exchanges.  Unfor¬ 
tunately,  there  are  no  ways  to  test  these  models  under  realistic  condi¬ 
tions,  i.e.,  an  atmospheric  nuclear  test.  Lven  the  use  of  high  explosive 
tests  has  been  criticized  as  not  simulating  a  nuclear  explosion  as  far  as 
the  dust  environment  is  concerned.  On  the  other  hand,  there  presently 
exists  a  large  data  base  from  the  atmospheric  test  program  dating  back  to 
the  1950s.  Much  of  these  data  were  recorded  on  photographic  films  and 
consist  of,  typically,  some  object  being  photographed  while  exposed  to 
the  effects  of  an  atmospheric  nuclear  burst.  For  many  of  the  films,  the 
data  include  the  shock-entrained  dust  (also  known  as  sweep-up  dust)  that 
was  generated  outside  of  the  crater  by  the  passage  over  the  ground  of  the 
shock  wave.  If  one  could  analyze  this  dust  versus  time  data  and  extract 
quantitative  dust  densities,  then  the  resulting  information  would  be 
very  useful  in  testing  different  computer  models  under  different  burst 
conditions . 

The  Contrast  Attenuation  Technique  (CAT)  for  deducing  dust  density 
from  photographic  film  was  developed  to  use  in  analyzing  existing 
films  (Reference  1).  Results  have  been  very  encouraging  and  a  considerable 
amount  of  data  have  been  and  continue  to  be  extracted  from  the  extensive  nuclear 
test  film  library.  The  purpose  of  the  MILL  RACE  validation  experiment  of 
the  CAT  methodology  was  to  reduce  the  uncertainties  in  using  these  data 
by  demonstrating  in  a  simple,  carefully  controlled  test,  the  sensitivity  of 
the  various  elements  used  in  the  CAT. 


1.  "Determination  of  Shock-Entrained  Dust  Concentration  from  Photographic 
Records  of  Nuclear  Weapon  Tests,"  R.  Snow,  M.L.  Price,  and  J.P.  Dotv, 
DNA  5119F  (Tune  1980). 
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A  theoretical  basis  of  CAT  is  given  in  Appendix  A;  as  derived, 
the  governing  equation  for  deducing  dust  de  ''ity  is  given  below: 


3 


__ 

z:  »2  f<“>  <■» 


o-i) 


where  is  the  dust  mass  concentiation ,  p  is  the  particle  grain  density, 
a  is  the  radius  of  the  dust  particle,  f(a)  is  the  particle  size  distribu¬ 
tion  function,  and  Q(a)  is  the  Mie  parameter.  The  remaining  symbol,  3, 
denotes  the  attenuation  coefficient  due  to  atmospheric  dust  associated 
with  the  contrast  between  a  target  and  its  background;  it  derives  from 
Koschmeider* s  Equation,  viz., 

C  =  Cq  exp(-BR),  (1-2) 

where  C,  C^  are  the  apparent  and  inherent  contrast,  respectively,  and  R 
is  the  distance  between  the  target  and  the  observer.  The  MILL  RACE  CAT 
test  consisted  of  six  spaced  photo  targets  at  both  the  nominal  10  psi  and 
20  psi  distance  from  GZ.  In  addition,  there  were  three  movie  cameras  at 
each  location  photographing  the  targets  from  just  before  detonation  until 

several  minutes  after.  By  using  Eqn.  (1-2)  for  several  different  targets, 
we  can  derive  the  following  basic  equation: 


til 


3(t)  = 


C.(t) 

cTTtY 


(1-3) 


R.  -  R. 
1  i 


where  i  and  j  denote  different  targets  and  the  range  Rj  is  larger  than  R^  . 
The  films  of  the  six  targets  were  analyzed  by  densitometi'y  and  the  contrast 
for  a  given  target  at  a  particular  time  deduced  by  using  the  following 
equation: 


C.(t) 


W, 


(1-4) 


where  and  H  refer  to  the  irradiance  of  the  film  due  to  the  background 
and  target,  respectively,  and  are  found  from  the  measured  density,  D,  of  the 
developed  film  by  use  of  the  calibrated  H  &  1)  curve  associated  with  the 
particular  film. 

Once  a  temporal  record  has  been  found  of  at  least  two  different 
photo  targets,  then  Eqn.  (1-3)  can  be  used  to  obtain  a  temporal  record  of 
6.  This  was  the  fundamental  objective  of  the  CAT  validation  experiment. 


Referring  to  Eqn .  (1-1),  the  determination  of  absolute  dust  densi¬ 
ties  requires  knowledge  of  the  particle  size  distribtuion  function  as  well 
as  the  Mie  parameter.  During  the  MILL  RACE  test,  a  separate  experiment 
fielded  by  SAI  (Reference  2)  measured  the  particle  size  distribution  vs 
time  and  the  absolute  dust  mass  concentration  vs  time  at  several  locations 
directly  in  the  line  of  sight  between  our  cameras  and  photo  targets. 

The  particle  size  distribution  function  was  supplied  to  us  by  SAI  and  we  calcu¬ 
lated  the  Mie  parameter  using  a  measured  index  of  refraction  value  for 
the  actual  dust  at  the  MILL  RACE  site.  The  grain  density,  P  ,  was  also 
measured  for  MILL  RACE  dust  and  along  with  our  deduced  value  of  3  completes 
the  terms  on  the  right  hand  side  of  Eqn.  (1-1);  this  immediately  gives 

p  (t).  We  can  then  compare  the  deduced  value  of  p  (t)  found  by  the  CAT 
m  m 

with  the  measured  value  of  p  ( t)  found  by  SAI  to  determine  the  accuracy 

m 

of  the  method. 

Section  2  of  this  report  describes  the  design  of  the  experiment 
while  Section  3  discusses  the  data  obtained  during  the  MILL  RACE  test. 

In  Section  4,  we  compare  the  dust  mass  density  deduced  from  the  CAT  with 
the  measured  data  supplied  by  SAI.  From  this  comparison  we  present  the 
conclusions  that  we  believe  are  warranted  from  the  experiment.  Several 
Appendices  are  attached  showing  the  theoretical  basis  of  the  technique 
in  more  detail. 
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"Measurements  of  Dust  Size  Distribution  and  Dust  Density  at  MILL  RACE," 
J.E.  Cockayne,  SAT,  (These  Proceedings). 
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Table  (2-1).  Definition  of  parameters  used  in  this  Section. 


PARAMETER 

T 

L 

S 

t 

F 

s 

f 

Y 

H 

D 

pps 


DEFINITION _ _ 

Film  running  time 
Film  length 
Camera  speed 
Exposure  time 
Frame  rate 
Sector 

Focal  length 

Slope  of  the  emulsion  characteristic 
curve 

Film  exposure 
Film  density 
Pulses  per  second 
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EXPERIMENT  DESIGN 

2.1  OBJECTIVES 

The  Contrast  Attenuation  Technique  validation  experiment  was 
planned  to  duplicate  as  closely  as  possible  the  test  hardware  used  for 
many  of  the  NTS  atmospheric  shots.  Because  the  NTS  atmospheric  test 
series  ended  over  20  years  ago,  obviously  it  would  be  difficult  to 
assemble  the  original  personnel  and  equipment  in  the  same  configuration 
as  used  at  that  time.  After  some  research  and  many  helpful  discussions 
with  members  of  the  Los  Alamos  EG  &  G  office,  we  finally  selected  the 
company  of  HSS,  Inc.  in  Bedford,  Mass.,  to  act  as  our  sub-contractor  in 
all  of  the  experimental  areas.  This  choice  was  based  on  HSS's  intimate 
knowledge  of  NTS  photographic  data  through  actual  experience.  All  of 
the  key  staff  members  at  HSS  who  helped  us  in  the  MILL  RACE  experiment 
had  been  involved  in  on-site  test  participation  at  NTS  in  the  photo¬ 
graphic  experimental  area.  Net  only  did  they  have  actual  experience 
dating  back  to  the  atmospheric  nuclear  tests,  but  HSS  also  retained  many 
of  the  same  cameras  that  were  used  at  those  tests.  In  addition,  HSS 
personnel  have  been  activd  in  many  other  recent  photographic  tests 
involving  H.E.,  rockets,  satellites,  etc.,  so  their  expertise  continues 
into  the  present.  Based  on  their  performance  at  MILL  RACE,  we  feel  very 
fortunate  in  choosing  HSS  to  assist  us. 

The  experimental  objectives  of  the  MILL  RACE  validation  test  of 
the  CAT  were  to  obtain  photographic  data  similar  to  NTS  data  showing  the 
effects  of  shock-entrained  dust  obscuring  known  targets  while  at  the 
same  time  an  absolute  measurement  of  the  dust  density  was  obtained.  Then 
by  comparing  the  dust  densities  deduced  from  the  photographic  data  with 
the  measured  data,  it  would  be  possible  to  determine  the  accuracy  of  the 
CAT.  Originally,  TBE  planned  to  conduct  both  the  photographic  and 
measurement  tasks  at  MILL  RACE.  Because  of  several  different  reasons, 
the  two  tasks  were  split  by  DNA;  SAI  in  McLean,  VA.  was  contracted  to 
perform  the  dust  measurements  while  TBE  conducted  the  photographic  portion. 

Data  from  previous  high  explosive  tests,  e.g..  Dice  Throw  and  Misers 
Bluff,  were  used  to  initially  specify  the  test  environment  for  MILL  RACE. 
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Based  on  movies  and  data,  we  selected  the  locations  for  our  two  test 
stations  at  the  nominal  10  and  20  psi  distances  from  GZ.  Originally,  we 
were  planning  to  use  three  test  stations — one  additional  at  15  psi — but 
finally  settled  on  two  to  minimize  the  number  of  dust  measurement  instru¬ 
ments  which  would  be  necessary.  Since  we  planned  to  use  six  cameras  for 
the  three  stations,  we  decided  to  continue  with  six  even  for  two  stations 


2.2  TEST  SITE  AND  EXPERIMENT  LOCATIONS 

The  MILL  RACE  shot  was  held  on  the  White  Sands  MissiLe  Range 
about  twenty  miles  south  of  the  Stallion  Range  Center  which  is  at  the 
northern  boundary  of  the  Base.  Starting  at  ground  zero,  the  terrain  had 
a  gentle  slope  rising  towards  the  east.  The  Contrast  Attenuation  Techni¬ 
que  Experiment  was  located  in  the  northeast  quadrant  approximately  between 
the  012  and  028°  radial.  Two  different  test  stations  were  placed  on  the 
test  bed — one  at  the  MO  psi  distance  (228.6m)  and  the  other  at  the  MO  ps  i 
distance  (298.7m).  Figure  (2-1)  shows  a  photo/layout  of  the  20  psi  test 
station;  the  10  psi  station  is  shown  in  F  ,ure  (2-2).  In  both  stations, 
the  cameras  were  aimed  towards  the  east 

A  test  station  consisted  of  two  camera  enclosures  and  six 
photo-targets.  Since  cameras  have  a  straight  line-of-sight .  we  obviously 
could  not  maintain  exactly  the  same  distance  from  GZ  for  all  of  the  units. 
Therefore,  we  chose  the  camera  enclosure  containing  the  GSAP  camera  and 
the  (1.83  x  1.83m)  photo-target  located  27.43m  from  the  camera  to  be  on 
the  20  psi  or  10  psi  distance.  Using  these  reference  points,  the  other 
five  photo-targets  were  located  as  shown  in  Table  (2-2). 

The  original  design  of  the  photo-targets  and  their  locations  was 
meant  to  provide  a  field-of-view  to  the  GSAP  cameras  as  shown  in  Figure  (2-3) . 
As  shown,  the  horizontal  spacing  between  targets  was  equal  to  the  target's 
angular  size  while  the  vertical  spacing  was  one-half  the  target's  angular 
size.  These  dimensions  were  modified  slightly  because  of  the  steel  columns 
supporting  the  plywood  sections.  There  was  no  impact  of  these  columns  on 
any  data  or  analysis. 

Height  of  the  field-of-view  was  M.52m  above  the  ground.  Since 
there  was  a  slight  upward  slope  to  the  ground  as  one  moved  from  the  cameras 
to  the  photo-targets,  this  meant  a  slight  elevation  to  the  field-of-view; 
the  actual  slope  of  the  field-of-view  was  about  two  degrees  for  both 
stations.  A  final  measurement  was  made  of  all  of  the  relevant  dimensions 
and  these  are  shown  in  Table  (2-3). 
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Fig. (2-1)  Photo  showing  layout  of  the  20  psi  test  station. 
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Fig. (2-2)  Photo  showing  layout  of  the  10  psi  test  station 


Table  (2-2) 


Location  of  photo-targets  on  the  MILL  RACE 
test  bed.  Also  shown  are  the  differences 
in  distance  from  GZ  and  the  estimated  dif¬ 
ference  in  overpressure  from  the  nominal 
10  or  20  psi 


Distance 

from 

camera 

(m) 

- - - - 

Size 

(m) 

Nominal  Range/ 
overpressure 

R/P  (m/psi) 

AR/AP 
(m/psi ) 

3.05 

.20  X  .20 

228.60/20 

+  .57/-. 15 

5.18 

.34  x  .34 

228.6/20 

-  .25/+. 07 

9.14 

.61  x  .61 

228.6/20 

-  1 .58/+. 41 

15.24 

1.00  x  1.00 

228.6/20 

+  2. 42/-. 62 

27.43 

1.83  x  1.83 

228.6/20 

0/0 

45.72 

3.05  x  3.05 

228.6/20 

-  4.21/+1 .07 

3.05 

.20  x  .20 

298.7/10 

+1.53/-. 05 

5.18 

•34  x  .34 

298.7/10 

-0 . 19/+.02 

9.14 

.61  x  .61 

298.7/10 

-1 .50/+. 15 

15.24 

1.00  x  1.00 

298.7/10 

+2. 32/-. 23 

27.43 

1.83  x  1.83 

298.7/10 

0/0 

45.72 

3.05  x  3.05 

298.7/10 

-4. 66/+. 45 

i 
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Table  (2-3)  Measured  distances  to  the  photo-targets  from  the  lens  of  the 
GSAP  camera. 


20  psi  Station 


Target 

Distance  from  GSAP 

Height 

above  ground 

(m) 

(m) 

1 

2.92 

1.45' 

2 

5.05 

1.39 

Top  Of 
target 

3 

9.20 

1.22 

4 

15.16 

1.58' 

5 

27.26 

1.68 

Bottom  of 
target 

6 

45.74 

1.65 

10  psi  Station 


7 

2.90 

1 .44 

8 

5.05 

1.29 

9 

9.14 

1.20 

10 

15.26 

1.60 

11 

27.00 

1.58 

Top  of 
target 
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45.17 


Bottom  of 
target 


By  using  the  GSAP  camera  location  as  the  reference  point  for 
locating  the  photo-targets,  the  f ield-of-view  for  the  other  two  cameras 
at  the  same  station  will  be  skewed  from  that  of  the  GSAP's.  For  the 
Flight  Research  camera,  the  distance  between  it  and  the  GASP's  was  about 
twenty  cm  so  that  their  f ields-of-view  were  only  slightly  different. 

The  Flight  Research  camera  was  able  to  see  all  of  the  photo-targets  with 
sufficient  spacing  between  them  for  contrast  measurement. 

However,  the  Mitchell  camera's  line-of-sight  was  displaced  later 
ally  almost  1.22  meters  from  the  GSAP's.  Due  to  this,  the  Mitchell's 
line-of-sight  had  almost  total  vignetting  of  the  2nd  and  3rd  targets  by 
the  1st  target;  the  farthest  three  targets  were  not  affected  to  any 
appreciable  degree. 
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2.3  PHOTO- TARGETS 

The  photo-targets  consisted  of  twelve  plywood  squares  (six  each 
at  the  20  psi  and  the  10  psi  stations)  and  of  dimensions  and  distance 
from  the  cameras  as  shown  in  Table  (2-2).  Although  both  stations  used 
identically  sized  wooden  squares  for  photo-targets,  the  frames  and 
supporting  columns  were  constructed  from  different  size  material  depend¬ 
ing  upon  whether  the  target  was  placed  at  the  20  psi  or  the  10  psi 
distance.  Because  damage  occurred  to  the  largest  target  (3.05m  x  3.05m) 
at  the  20  psi  station,  a  description  of  that  particular  unit  will  be 
given  b.  low. 

The  frame  for  the  large  plywood  section  was  made  from  12.7cm  (5") 
nominal  ID  standard  weight  steel  pipe  with  5.08  x  5.08cm  (2"x  2")  x  .48cm 
(3/16")  angle  legs  welded  to  the  pipe  in  sections  of  from  30cm  to  68cm  in 
length.  Prior  to  the  detonation,  an  examination  of  this  photo-target  (as 
well  as  the  others)  indicated  very  close  alignment  of  the  radial  emanating 
from  ground  zero  with  the  plane  of  the  plywood  target.  If  a  perfectly 
normal  shock  wave  had  arrived  at  this  target,  the  only  cross-section  to 
intercept  the  pressure  wave  consisted  of  the  5"  ID  pipe  on  the  windward 
column;  therefore,  one  would  expect  any  damage  to  be  in  a  radial  direction 
away  from  GZ. 

All  of  the  photo-targets  and  supporting  members  were  painted  prior 
to  the  explosion  with  a  high-quality  flat  black  paint.  The  intent  was  to 
make  the  targets  appear  to  be  ideal  black  bodies  throughout  the  spectral 
range  of  the  film.  Figure  (2-3)  shows  the  photo-targets  in  the  field-of- 
view  of  the  GSAP  cameras  prior  to  the  test.  Unfortunately,  events  occur¬ 
ring  during  the  actual  HE  test  severely  modified  our  carefully  prepared 
surface  and  eliminated  the  black  body  hypothesis. 


2.4  CAMERAS 

Three  different  camera  tvpes  were  used  at  both  the  20  psi  and  the 
10  psi  station.  A  Mitchell  HS-100  high-speed  35mm  camera  and  a  Kodak  GSAP  (Oun 
Sight  Aiming  Point)  16mm  camera  were  chosen  as  identical  to  cameras  used 
during  NTS  shots.  In  addition,  a  Flight  Research  IV  C  35mm  camera  was  also 
used  for  several  reasons: 

1.  To  provide  redundancy  for  the  Mitchell  camera. 

2.  To  obtain  data  on  a  more  recent  vintage  camera. 

3.  To  obtain  35mm  data  at  approximately  the  same  location  as 
the  16mm  GSAP  in  order  to  compare  film  characteristics. 

All  cameras  were  set  with  a  lens  focus  at  infinity  for  all  runs. 

2.4.1  Mitchell  35mm  Camera 

The  Mitchell  cameras,  made  twenty-five  or  thirty  years  ago,  are  large, 
high  quality,  professional  movie  cameras  that  have  complex  mechanisms  and  are 
powered  by  110V  A.C.  electric  motors.  Film  magazines  for  this  camera  handle 
up  to  152m.  film  lengths  which  was  longer  than  the  other  two  camera's  maximum 
film  length.  Film  running  time,  T,  can  be  calculated  as  follows: 


where  L  is  the  film  length  in  meters,  S  is  the  camera  speed  in  frames/sec, 
and  F  is  the  number  of  frames  per  meter.  For  35mm  film  F  =  52.5.  Thus 
for  the  Mitchell  cameras, 

_  8  x  103  , 

T  =  — — —  seconds, 

or 

o,  frames 

T  =  333  seconds  at  S  =  24  — — - , 

sec 

the  speed  used  during  the  MILT.  RACE  tests.  Assuming  a  nominal  value  of  10 
seconds  to  obtain  the  24  frames/sec  speed  after  turning  the  Mitchell's 
motor  on,  we  find 

T  =  323  seconds  =  5  min  .’3  sec 

as  the  total  film  time  available  for  data  gathering.  (Actually  the  test 
time  was  slightly  greater  than  this  because  the  camera  was  operating 
at  less  than  24  frames/sec  during  the  initial  10  seconds  i.e.,  accelerating 
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from  zero  to  the  final  speed.  Therefore  the  amount  of  film  exposed 
during  this  10  seconds  was  less  than  the  5  meters  calculated  using  a 
constant  speed  of  24  frames/sec  and  thus  there  was  more  film  available 
during  the  data  gathering  time) . 

An  important  feature  on  the  Mitchell  cameras  not  included  on 
either  the  Flight  Research  or  the  GSAP  cameras  were  timing  marks  placed 
on  the  film  at  a  50pps  rate.  A  TD-60  Timing  Mark  Generator  was  used  to 
provide  accurate  timing  signals  which  were  then  used  to  fire  spark  markers. 

The  timing  mark  was  located  on  the  outer  edge  of  the  Mitchel  film  and 
after  processing  provided  a  means  for  determining  the  time  at  which  an 
event  occurred  on  the  film.  The  use  of  spark  markers  instead  of  neon 
bulbs  was  mandated  by  the  low  ASA  number  of  the  3414  film.  During  pre¬ 
liminary  tests  at  HSS ,  Inc.,  the  neon  markers  normally  used  on  the  Mitchell 
cameras,  were  not  able  to  expose  the  film  sufficiently  and  thus  were  not 
visible  after  processing.  Not  coincidentally  ,  during  the  atmospheric  nuclear 
testing  t  an  0f  the  Mitchell  cameras  were  adapted  to  spark  marker  systems 
in  order  to  expose  the  MF  film  used  for  most  of  the  shots  because  it  also  had 
a  very  low  ASA  number. 

Exposure  of  film  inside  a  Mitchell  camera  can  be  controlled  using 
several  independent  settings.  The  frame  rate,  F,  shutter  opening  (f/stop), 
and  sector,  s,  all  determine  the  exposure.  The  exposure  time,  t,  is 
given  by 


t  =  -(-?—) 

FV360; 


(2-2) 


For  the  MILL  RACE  shot,  the  frame  rate  was  24  frames  and  the  sector  was 


170  .  Therefore 

t  n  . 02  sec  =  ^  sec . 

Both  Mitchell  cameras  used  a  35mm  Baltar  lens. 


sec 


The  f/stop  determined  by  pre-shot  testing  was  f/3.5  for  the  10  psi  Mitchell 
and  f/2.8  for  the  20  psi  Mitchell. 

2.4.2.  Kodak  GSAP  16mm  Cameras 

The  Gun  Sight  Aiming  Point  cameras,  manufactured  by  Kodak,  were 
of  WW1I  vintage.  Due  to  their  small  size,  ruggedness,  low-cost,  and 
simplicity,  a  large  number  were  used  during  the  atmospheric  nuclear  test¬ 
ing  program.  Most  of  the  NTS  film  library  consists  of  16mm  records  taken 
with  GSAP  cameras.  Therefore,  we  placed  great  emphasis  on  obtaining  data 


using  GSAP  cameras  at  MILL  RACE.  The  cameras  themselves  were  veterans 
of  several  NTS  shots.  As  expected,  these  cameras  were  trouble-free  and 
produced  relatively  high  quality  data  in  spite  of  their  age.  Because  we 
were  unable  to  procure  the  Type  3414  Film  in  16mm  size,  we  used  Kodak 
Type  2476  film  in  15.2m  rolls  which  had  the  same  spectral  response  but 
a  much  higher  ASA  number.  Both  of  the  GSAP  cameras  were  run  at  16  frames/sec 


with  a  sector  of  130 


The  corresponding  exposure  time,  t,  is  given  by 


1,  s  1  , 130.  _ 

C  F(360)  16  360  *°2  Sec> 


or 


t  =  sec;  the  same  as  the  Mitchells’. 

Since  16mm  film  has  131.2  frames/m,  the  total  running  time  for  the  GSAP 
cameras  was 

T  =  ~  =  (1^ >2)  (131.2) _  _  125  seconds. 

S  16 

Unlike  the  Mitchell  cameras,  the  GSAP  required  a  very  short  acceleration 
time  so  that  the  available  data  gathering  time  was  practically  identical 
to  T.  There  were  no  timing  marks  placed  on  the  GSAP's  film  during  the 
run  but  speed  tests  conducted  prior  to  MILL  RACE  determined  the  fixed 
framing  rate  for  both  cameras.  Therefore,  it  was  simply  a  matter  of 
measuring  the  distance  along  the  film  to  determine  the  exposure  time  of 
any  particular  event.  The  lens  used  in  the  GSAP  cameras  were  Bell  and  Howell 
Super  Comet  with  a  focal  length  of  17.8mm.  During  MILT-  RACE  the  GSAP's  were 
run  using  f/stops  of  f/19  and  f / 1 6  for  the  10  psi  and  20  psi  cameras,  respec¬ 
tively. 


2.4.3.  Flight  Research  Cameras 

The  Flight  Research  35mm  cameras  are  of  fairly  recent  design,  and 
classify  as  relatively  simple  units  about  one-third  the  size  of  the  Mitchell 
camera.  Because  of  their  compact  size,  it  was  possible  to  place  the  Flight 
Research  cameras  in  the  same  enclosure  as  the  GSAP  cameras  so  that  their 
f ields-of-view  were  essentially  identical.  This  allowed  for  a  direct  com¬ 
parison  during  film  analysis  of  the  Type  3414  and  2476  Kodak  films. 

Framing  of  the  camera  was  controlled  by  a  fixed  speed  motor  and 
gear  selection  for  speed  control.  During  the  MILT-  RACE  tests,  they  were 
run  at  10  frames/sec  and  set  for  a  sector  of  72°.  Thus  the  exposure  time. 
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t,  is  given  by 


1  /  s  \  X  /  7  2  »  « « 

C  =  T"  (360)  =  ~W  (  360  "  -°2 


or 


t  = 


50 


sec. 


The  data  gathering  time,  T,  is  found  from 


T  =  LF 


shot  runs  at  MILL  RACE  used  21.3m  rolls  of  3414  film  so 
(21.3)  (52.5) 


T  = 


10 


=  112  sec. 


The  Flight  Research  cameras  used  35mm  Camera  Equipment  Corp.  lens  with 
an  f/stop  of  f/3.5  for  both  the  10  psi  and  20  psi  station. 


2.5  FILM 

The  film  used  during  the  NTS  nuclear  testing  was  Kodak  Type  Micro¬ 
file  for  most  of  the  shots.  This  film  was  especially  specified  for  use  in 
a  nuclear  environment  because  of  its  low  sensitivity  to  gamma  radiation. 

One  side  effect  of  such  a  specification  was  a  very  low  ASA  rating  for 
this  film  (ASA  1-3)  .  Kodak  no  longer  manufacturers  this  particular 
emulsion  so  an  extensive  investigation  was  made  of  the  available  films 
to  find  the  best  substitution.  Aiding  in  this  search  was  Mr.  Charles  Wyckoff 
who  was  a  major  contributor  to  the  initial  specification  of  the  MF  film.  He 
recommended  Kodak  Type  3414  film  as  the  best  substitute  for  the  MF  film  based 
on  very  similar,  if  not  identical,  spectral  response  functions.  The  major 
differences  between  the  two  films  was  apparent  in  the  curves  showing  film 
density  vs  exposure  for  various  choices  of  y  (the  slope  of  the  emulsion  char¬ 
acteristic  H&D  curve);  the  MF  film  had  a  much  larger  latitude.  Mr.  Wyckoff 
believed  that  this  difference  would  be  minimized  by  using  suitable  processing 
techniques  and  would  not  pose  any  problem  in  analyzing  the  film  data  for  our 
tests . 

Data  for  type  3414  film  is  shown  on  the  accompanying  data  sheets 
supplied  by  Kodak  (see  Appendix  C) .  Upon  selection  of  this  film,  we  con¬ 
tacted  Kodak  for  price  and  delivery  information.  The  result  of  our  inquiry 
was  an  order  for  sixteen  rolls  of  35mm  film  in  152.4m  lengths.  It  was 

not  feasible  to  purchase  any  16mm  film  with  3414  emulsion  and  for  a  while 

we  were  even  apprehensive  about  obtaining  the  35mm  film  in  time  for  MILL 
RACE.  This  last  fear  proved  negative  and  we  received  the  ordered  film  a 
month  prior  to  the  MILL  RACE  fielding. 

However,  we  still  had  to  choose  the  16mm  film  for  use  in  the  GSAP 

cameras.  Again  Mr.  Wyckoff  assisted  in  the  selection  and  recommended  use 

of  Kodak  Shellburst  Film  #2476.  As  evidenced  on  the  accompanying  data 
sheets,  the  main  difference  between  3414  and  2476  lies  in  the  exposure 
index  with  values  of  2.5  and  250,  respectively.  Their  spectral  response 
functions  although  not  identical,  are  very  similar. 

After  deciding  upon  2476  film  for  our  use,  we  contacted  Kodak  only 
to  learn  that  they  were  completely  out-of-stock  for  this  film  and  were  not 
planning  to  manufacture  additional  units  for  six  months.  With  the  assis¬ 
tance  of  the  MI’LL  RACE  Test  Staff,  we  finally  located  and  procured  sufficient 
quantities  of  16mm  2476  film  through  the  WSMR  photographic  department. 
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POWER  AND  CONTROL 

The  power  and  control  system  used  for  the  CAT  validation  experiment 
is  shown  in  Fig.  (2-4).  Except  for  the  two  camera  enclosures,  the  re¬ 
maining  units  were  placed  inside  an  excavated  pit  next  to  the  enclosures 
and  covered  with  sandbags  prior  to  detonation.  A  number  of  foam  sheets 
were  placed  under  and  around  the  units  to  provide  some  shock  protection. 

Power  to  the  cameras  was  provided  by  commercial-type  automobile 
batteries  newly  procured  for  the  test.  Battery  charging  was  performed 
as  necessary  and  was  accomplished  by  having  two  sets  of  batteries,  one 
being  charged  while  the  other  was  used. 

The  cameras  were  turned  on  by  a  relay  closure  circuit  provided 
by  Field  Command.  This  signal  was  used  to  activate  our  control  unit  which 
had  a  latching  relay  to  keep  power  applied  to  the  cameras  even  after  the 
relay  closure  provided  by  Field  Command  was  disabled.  Thus  we  planned  to 
activate  our  cameras  beginning  at  10  seconds  prior  to  detonation  and?$Leep 
them  on  until  we  arrived  at  the  test  stations  and  manually  disabled  them. 

The  Mitchell  cameras  had  a  built  in  timing  mark  generator  to 
provide  a  known  reference.  This  unit  was  controlled  by  a  TD-60A  film 
marker  set  to  operate  at  50  pps. 
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SECTION  3 
DISCUSSION 

3.1  TEST  RESULTS 

The  MILL  RACE  explosion  occurred  at  about  1235  on  Sept.  16,  1981. 

All  of  our  cameras  operated  as  planned  except  for  the  Mitchell  at  10  psi 
which  jammed  after  about  one  half  of  the  film  was  exposed.  Damage  occurred 
to  the  largest  photo-target  at  the  20  psi  station;  also,  the  Mitchell  cam¬ 
era  enclosure  at  20  psi  was  slightly  tilted  due  to  the  force  of  the  shock 
wave,  but  a  boresight  check  of  the  camera  indicated  that  all  of  the  photo¬ 
targets  were  still  within  its  f ield-of-view.  Post-test  examination  revealed 
that  all  the  photo-targets  at  the  20  psi  station  and  the  three  smallest  at 
the  10  psi  station  were  coated  with  a  thin  layer  of  dust  which  was  embedded 
on  the  painted  surface.  It  was  obvious  that  these  targets  no  longer  approx¬ 
imated  an  ideal  black  surface.  By  using  two  calibrated  reflection  cards , 
we  estimated  the  surfaces  to  have  a  reflection  of  between  (0.2-0. 3). 

Fig.  (3-1)  and  (3-2)  show  post-test  photos  of  the  photo-targets  at  the  two 
test  stations. 

During  removal  of  the  equipment  after  the  experiment  was  completed, 
we  discovered  the  TD-60A  units  at  the  20  psi  station  had  received  a  sizable 
shock  which  resulted  in  the  disconnection  of  several  circuit  boards.  Most 
likely,  the  ground  shock  from  the  detonation  caused  this  and  ended  the 
timing  mark  generation  at  about  shock  arrival  time. 

The  data  film  were  caiefully  removed  and  packed  for  shipment  to 
the  processing  laboratory.  Under  the  supervision  of  Mr.  Charles  Wyckoff, 
a  step  wedge  filter  was  used  to  expose  identical  film  and  along  with  the 
data  film,  processed  according  to  the  technique  indicated  below. 
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3.2  FILM  PROCESSING 

3.2.1  Sensitometry  and  Film  Development 

3. 2. 1.1  Processing  (General) 

Eight  shot  films  of  two  different  black  &  white  film  stocks  were 
processed  by  Film  Facts  Engineering  (FFE)  Inc.  in  a  Kodak  Versamat  11C-M. 
All  films  were  proceeded  by  and  accompanied  by  sensitometric  step  control 
scales  using  identical  emulsion  stock.  An  EG&G  Xenon  Flash  Sensitometer 
was  used  to  expose  the  control  and  check  scales. 

The  Kodak  Versamat  11C-M  is  a  continuous  processing  eleven  inch 
wide  film  machine  with  a  large  variety  of  processing  controls  such  as 
temperature,  speed  drive  variation,  continuous  heating  and  cooling  devel¬ 
oper  temperature,  balancing  of  the  developer  with  constant  filtering,  opera 
tor's  option  of  replenishing  the  developer  and  fixer  solutions  with  fresh 
solutions  before,  during  or  after  processing  in  accordance  with  the  amount 
of  film  processed. 


3.2. 1.2.  Sensitometr 


FFE,  Inc.  used  an  EG&G  Xenon  Flash  Sensitometer  Mark  VII  to  expose 
a  calibrated  step  scale  on  to  fresh  film  stock  of  the  same  emulsion  batch 
used  for  the  record  film.  A  test  scale  was  processed  before  any  record 
film  was  run  or  processed  in  the  machine.  The  test  scale  was  read  on  a 
Macbeth  Quantalog  Densitometer.  Measurements  were  plotted  to  determine 
if  the  processing  machine  controls,  i.e.  ,  speed,  temperature,  drying, 
fixation  and  washing  were  all  satisfactory  for  processing  the  film.  When 
the  record  film  was  to  be  processed,  a  scale  was  placed  at  the  head  and 
at  the  tail  of  the  film.  It  was  rarely  possible  to  place  the  control 
scale  directly  on  a  shot  film.  An  alternative  was  generally  taken;  con¬ 
trol  scales  were  made  and  run  at  the  head  and  tail  of  the  record  film. 

Two  records  were  kept  of  the  machine  processing  conditions;  one  in  a  pro¬ 
cessing  1 og  in  the  machine  room  and  the  other  on  the  plot  of  the  control 
scale  (density  vs  log-exposure).  The  results  of  measurements  taken  for 
the  two  Flight  Research  and  one  GSAP  cameras  are  shown  in  Figs.  (3-3), 
(3-4)  and  (3-5).  (Note  that  the  units  shown  along  the  abscissa  are  not 
absolute  exposure  units  but  merely  abitrary  units  that  were  useful  for 
film  calibration.) 


3 . 2 . 1 . 3  Film  Development 

There  were  two  types  of  shot  films  used:  Kodak  High  Definition 
Aerial  Film  (3414-1751021  35mm),  Estar  base;  and  Kodak  Linagraph  Shellburst 
Film  (2476-49335  E  16mm)  Estar  base.  Kodak's  published  data  for  processing 
of  these  films  in  a  Versamat  processor  was  limited  to  very  specific  proces¬ 
sing  conditions.  For  example:  the  3414  emulsion  characteristic  curves  were 
specific  to  Versamat  885  chemicals,  one  developer  rack  and  speeds  of  10  ft/min 
and  15  ft/min;  no  emulsion  characteristic  curves  were  given  for  the  type  2476 
using  Versamat  processing;  however,  the  gamma  of  the  curves  was  given  for 
three  Versamat  chemistries  Type  A,  Type  885  and  Type  641.  The  information 
supplied  by  Kodak  was  of  little  value  because  their  processing  specifications 
led  in  every  case  to  high  contrasts  and  short  exposure  latitudes  for  both 
films.  As  a  result,  it  was  necessary,  under  the  guidance  of  Mr.  Charles  W. 
Wyckoff — a  recognized  expert  in  film  processing — to  experiment  with  various 
processing  chemistries  and  processing  conditions. 

The  processing  machine  options  available  required  both  close  moni¬ 
toring  and  regular  sensitometric  checks.  Many  different  developers  and 
processing  conditions  were  tested  to  achieve  optimum  performance  for  a  given 
application  with  a  unique  combination  of  developer,  time  and  temperature. 
Variance  from  these  conditions  will  yield  less  than  optimum  results  due  to: 
alteration  of  base  fog  level,  loss  of  film  sensitivity,  increase  in  grain 
size,  incomplete  development,  chemical  residue  left  in  the  film,  or  improper 
drying. 

The  optimum  chemistry  to  achieve  the  desired  MILL  RACE  Film  char¬ 
acteristics  was  found  to  be  Kodak  Type-B  Chemistry  with  a  slightly  modified 
pH  value.  Two  developer  racks  were  required.  Optimum  developing  temperature 
was  87°F.  Optimum  film  speeds  differed  for  the  two  emulsions.  Exact  proces¬ 
sing  conditions  for  each  of  the  eight  films  are  given  in  Table  (3-1). 

Adjacency  and  Bromide  Drag  effects  are  problems  that  constantly 
trouble  photometric  photography.  The  adjacency  effect  is  a  problem  associ¬ 
ated  with  two  or  more  images  which  lie  close  together.  Prolonged  or  irregu¬ 
lar  development  can  cause  loss  in  density  of  one  or  both  images  or  a  buildup 
of  density  between  the  two  images.  This  effect  was  minimized  or  eliminated 
in  the  Versamat  11C-M  by  proper  maintenance  of  the  machine  and  development 
of  the  film  in  the  mid-ranges  of  the  drive  speed  which  goes  from  0  to  25 


feet  per  minute.  Tri  addition,  the  machine  had  certain  of  the  rollers 
slightly  overdriven  to  eliminate  film  slack  due  to  normal  film  stretching. 
Bromide  buildup  in  the  developer  is  caused  by  an  imbalance  of  excessive 
amounts  of  image  processing  in  relation  to  amount  of  replenisher  being 
fed  to  the  developing  tank.  Bromide  drag  problems  were  eliminated  in 
this  machine  through  proper  replenishment  techniques,  recirculation  of 
the  developer,  filtering  of  the  developer,  and  sensitometric  checks  on 
development. 

Agitation  of  the  Versamat  11C-M  was  derived  from  the  movement  of 
the  film  itself  through  the  developer  tank,  assisted  by  chain-driven  rollers 
which  in  turn  drove  or  propelled  it  and/or  churned  the  adjacent  developer 
to  drive  fresh  developer  on  to  the  surface  and  at  the  same  time  drive  the 
used  developer  away  from  the  film  emulsion  surface.  The  only  other  motions 
in  the  developer  tank  were  the  recycling  of  the  developer  exiting  to  the 
cooling  and  heating  chambers  to  maintain  constant  temperature  of  the  eight 
gallons  in  the  machine  processing  tank,  and  a  trickle  input  of  fresh 
replenisher  developer  stock.  The  machine  operator  had  mechanical  controls 
of  this  vital  process  for  heating  and  cooling  of  the  developer  and  to  control 
the  drive  speed  of  the  machine  for  the  time  of  film  development.  These 
controls  yielded  the  total  development  gamma  resultant  on  the  film. 


3.2.2  Film  Densitometry 

After  all  of  the  films  were  processed,  we  spent  considerable  time 
examining  each  film  using  a  Moviola  film  display  unit.  All  three  cameras 
at  each  station  had  good  data  showing  the  effects  of  dust  on  the  contrast 
of  the  photo-targets.  Because  the  Mitchell  camera’s  f ield-of-view  did  not 
clearly  delineate  the  three  closest  targets  and  their  backgrounds,  and  due 
to  the  loss  of  timing  marks  at  shock  arrival  time,  we  decided  to  use  the 
Flight  Research  35mm  films  as  our  prime  source  of  data  for  densitometry 
measurements.  The  GSAP  cameras  would  be  used  at  selected  points  to  obtain 
data  on  the  16mm  2476  film  for  correlation  with  the  35mm  3414  film. 

Once  the  decision  to  use  the  35mm  Flight  Research  film  for  densito¬ 
metry  was  made,  we  chose  two  different  densitometers  to  use  for  the  analysis 
The  amount  of  data  available  was  considerable;  one  Flight  Research  camera 
alone  had  over  1000  frames  of  exposed  film.  Clearly,  we  needed  a  fairly 
quick  method  for  measuring  the  density  of  the  film  for  each  target  and  its 
associated  background.  However,  there  was  interest  in  more  detailed  study 
of  selected  frames  on  the  uniformity  of  the  dust  and/or  to  measure  the  dust 
gradient  across  an  individual  target.  To  accomplish  both  the  quick  analyses 
as  well  as  the  detailed  dust  profile  studies,  we  selected  two  different  den¬ 
sitometers:  the  Macbeth  Quantalog  Densitometer  provided  a  diffuse  density 

measurement  in  that  an  average  film  density  over  ~  1mm2  was  measured  using 
a  manually  controlled  scanning  system  for  both  the  target  and  its  background 
the  Photometric  Data  Systems  Model  -  1000  provided  a  microdensitometer  scan 
using  a  beam  with  an  area  of  =  4  x  10  3  mm2  and  a  digitally  controlled  scan¬ 
ning  motor.  More  information  is  provided  below  on  each  of  these  instruments 
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3 . 2 . 2 . 1  Diffuse  Density  Measurements  of  MILL  RACE  FILMS 

Diffuse  density  measurements  of  the  targets  and  backgrounds  on  the 
MILL  RACE  photographic  films  were  made  with  a  Model  TD-102  Macbeth 
Quantalog  Densitometer. 

All  diffuse  density  measurements  using  the  TD-102  were  made  with 
a  circular  aperture  0.046  inch  (1.17mm)  in  diameter.  This  aperture  was 
the  smallest  aperture  which  could  be  employed  with  the  TD-102  while  main¬ 
taining  a  valid  calibration  of  the  instrument.  Target  images  on  the  Flight 
Research  IV  C  and  Mitchell  HS-100  Camera  35mm  records  were  approximately  2  mm  x 
2  mm  in  size;  the  measuring  aperture  of  the  TD-102  could  comfortably  fit 
within  these  target  images.  On  the  other  hand,  target  images  on  the  GSAP 
Camera  records  were  about  1  mm  x  1  mm  in  size.  Because  the  GSAP  images 
were  smaller  than  the  densitometer  aperture  it  was  not  possible  to  make 
diffuse  density  measurements  on  the  GSAP  films  with  the  TD-102.  (Note: 
the  primary  reason  for  using  the  TD-102  rather  than  the  micro-densitometer 
for  the  contrast  measurements  was  that  the  TD-102  permitted  large  quantities 
of  data  to  be  processed  in  a  reasonable  amount  of  time.) 

Figure  (3-6)  illustrates  the  locations  on  the  Flight  Research  Cam¬ 
era  records  where  diffuse  density  measurements  were  made.  Background 
measurements  made  at  positions  Bl,  B2  and  B3  were  intended  for  use  with 
targets  T1 ,  T2  andT3,  respectively.  They  could  also  be  used  for  deter¬ 
mining  target  contrasts  with  targets  T4,  T5  and  T6  if  the  situation  permitted. 
At  late  times,  if  the  optical  thickness  of  the  atmospheric  dust  along 
lines-of-sight  Bl ,  B2  and  B3  was  not  infinite  (i.e.  desert  terrain  was 
showing  at  these  locations)  then  the  background  measurements  at  positions 
B4,  B5  and  B6  could  be  substituted  for  use  with  targets  T4 ,  T5  and  T6 
respectively  . 

Apparent  target  contrasts,  C,were  determined  in  the  manner  outlined 
in  Figure  (3-7).  An  emulsion  characteristic  curve  (H&D)  generated  from  a 
sensitometric  step  wedge  using  the  TD-102  densitometer  was  used  to  convert 
the  target  and  background  film  densities  to  relative  brightness  values. 

(The  sensitometric  step  wedge  was  processed  under  carefully  controlled 
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Figure  (3-7) .  Outline  of  procedure  for  obtaining  optical  extinction 
coefficients  of  atmospheric  dust  from  MILL  RACE  films. 
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conditions  along  with  the  data  film,  using  development  and  agitation 
techniques  designed  to  overcome  all  edge-effects  and  bromide-drag  effects). 

The  relative  brightness  values  were  then  used  to  determine  the 
apparent  contrast  of  any  given  targeL  using  as  the  defining  relation  for 


apparent  contrast 


where  is  the  relative  brightness  of  the  target  and  is  the  relative 
brightness  of  the  adjacent  dust  background. 

The  optical  extinction  coefficient  3  of  the  atmospheric  dust 
between  any  part  of  targets  was  determined  from  the  relation 


(1-3) 


where  i  designates  any  target  1  through  5,  and  j  designates  any  target  2 
through  6. 


3. 2. 2. 2  Microdensitometer  Scanning  of  MILL  RACE  Films 

To  supplement  the  diffuse  density  measurements  made  on  the  MILL 
RACE  films  with  the  Macbeth  Densitometer  and  to  examine  spatial  variations 
in  the  target  scenes,  scans  of  selected  frames  were  made  on  a  recording 
microdensitometer.  This  type  of  densitometer  measures  specular  density 
rather  than  diffuse  density  and  a  separately  determined  H&D  curve  mea¬ 
sured  by  the  instrument  itself  must  be  obtained.  While  the  corresponding 
density  measurements  between  the  two  methods  cannot  be  compared,  the 
relative  exposure  values  obtained  using  the  respective  H&D  curves  can  be 
directly  compared. 

The  microdensitometer  used  for  these  scans  was  Photometric  Data 
Systems  Model-1000  (PDS-1000)  .  In  this  instrument  the  film  was  mounted 
between  clear  glass  plates  which  sit  on  a  platen  driven  in  an  x-y  rectan¬ 
gular  coordinate  system  by  two  digital  step  motors.  The  minimum  step  in 
either  direction  was  2  microns.  The  motors  could  be  driven  continuously  or 
with  single  steps.  They  could  be  manually  operated  from  front  panel  switches 
or  automatically  programmed  by  an  auxiliary  control  system. 

A  variety  of  scanning  apertures  was  available  on  a  rotating  wheel 
A  71  micron  diameter  circular  aperture  was  used  in  all  of  the  scanning 
done  on  these  films.  The  optical  system  was  in  the  form  of  a  single  beam 
transmissometer  whose  output  voltage  was  converted  to  density  by  an  elec¬ 
tronic  logarithmic  amplifier.  The  resulting  analog  voltage  output  could  be 
recorded  on  a  chart  recorder  or  converted  to  digital  output  and  recorded 
on  a  magnetic  tape. 

For  these  scans,  the  analog  density  output  was  recorded  on  a  Brush 
Recorder  Mark  250  which  gave  an  ink  record  on  a  11.43cm  (4V)  wide  paper 
chart  with  an  imprinted  precision  rectilinear  grid  graduated  in  50  divis¬ 
ions  for  the  44"  width. 

The  experimental  procedure  started  by  mounting  the  film  on  the 
platen  between  clear  glass  plates  and  aligning  the  film  with  the  x-y  scan¬ 
ning  coordinate  system.  The  optical  system  was  focused  and  the  transmis¬ 
sion  output  was  set  to  read  unity  on  a  front  panel  digital  voltmeter  when 
the  light  beam  was  passing  through  clear  glass.  The  output  was  then 


switched  to  read  density  and  the  density  scale  was  set  using  c.i !  il>  c.U  <  d 
Density  filters  mounted  on  an  internal  rotating  filter  wheel  which,  inter 
cepts  the  optical  beams.  Next,  the  recorder  scale  was  adjusred  to  re-ui 
density  as  displayed  on  the  front  panel  digital  meter.  All  of  rhc 
electronic  settings  were  reset  before  each  scan.  The  platen  srann 
speed  and  the  chart  speed  were  chosen  to  give  a  suitable  record  longm 
The  recorder  amplifier  was  set  to  give  near  full  scale  density  readi'- -s 
for  each  scan.  When  the  density  varied  widely  during  a  scan,  tne  measn. 
ment  was  repeated  using  a  different  recorder  setting  Ln  ocdc>-  to  "h 
good  readings  at  all  points. 
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3.3  CONTRAST  MEASUREMENTS 

The  measured  film  density  found  usirg  the  densitometer  converts 
to  an  exposure  value  by  using  the  H  and  D  calibration  curve.  For  a 
given  film  frame,  each  target  that  was  visible  had  a  pair  of  exposure 
values  assigned  to  it — one  for  the  target  and  one  for  the  background. 
These  exposures  related  directly  to  the  apparent  brightness  of  the  target 
and  its  background.  Thus,  we  define  the  contrast  of  a  target,  T, 
at  distance  R  by  the  following  expression: 


C(R) 


(3-1) 


where  and  refer  to  the  brightness  of  the  background  and  the  bright¬ 
ness  of  the  target,  respectively.  Since  the  brightness  and  the  exposure 
were  linear  in  this  application,  we  substituted  one  for  the  other  in 
determining  the  contrast,  i.e., 

(3-2) 

The  contrast  for  a  given  target  at  a  given  time  was  found  by  using 
the  measured  exposure  value  for  the  target  and  its  background  along  with 
Eqn.(3-2).  Both  of  the  Flight  Research  films  were  analyzed  out  to  60 
seconds  after  shock  arrival.  An  analysis  shown  in  the  next  section  derives 
a  few  simple  selection  rules  used  for  deciding  which  targets  would  pro¬ 
duce  the  best  data  when  multiple  targets  were  visible. 

The  background  exposure  on  the  film  was  determined  by  many  factors. 
All  sources  of  light  such  as  direct  sunlight,  scattered  light,  reflected 
light,  and  emitted  light  from  any  object  contributed  to  the  back¬ 
ground  brightness.  An  equation  for  the  background  brightness,  B^ ,  is 
given  in  Eqn.(3-3). 

00 

B  =  /  H<t>oexp  fr(3o  r)]  dr. 


H,  -  H„ 


C(R)  = 


«b 


(3-3) 


where  H  is  the  total  airlight  at  a  distance  r,  (see  Fig.  (A-l),  i.  the 

total  scattering  coefficient,  is  the  angular  dependence  function  of 

and  8  is  the  extinction  coefficient.  When  H,  <p  and  8  are  uniform 
ext  ext 

between  the  camera  and  r,  we  can  integrate  Fqn.  (3-3)  and  obtain  the  fo 
lowing  term: 


Bb  = 


ext 


(3-4) 


Eqn.  (3-4)  applies  for  optically  dense  media,  suer  . s  heavy  dust  concen¬ 
trations,  when  (g  r)  >>1.  Even  when  Eqn.  (3-4)  can  be  used,  however, 
the  background  radiance  in  general  will  be  changing  witn  time  n«  H,  <p  and 
g  change. 

Contrast  was  defined  by  Eqn.  (3-1).  If  the  target  n r  ■  <  s  was 

zero,  which  would  occur  for  an  ideal  black  target,  then  the  in'-  -c.it  con¬ 
trast  would  equal! 


C 

o 


1* 


(3-5) 


for  all  values  of  Bb>  This  was  the  initial  objective  of  our  MILL  RACE 
experiment  design;  we  carefully  painted  all  of  our  targets  with  a  high 
quality,  flat  black  paint  in  order  to  satisfy  Eqn.  (3—5).  Unfortunately, 
the  surfaces  of  the  targets  were  significantly  modified  by  a  layer  of  dust 
embedded  into  the  paint  sometime  during  the  passage  of  the  dust  cloud. 

The  resulting  surface  no  longer  approximated  an  ideal  black  surface  but 
probably  had  an  refectivity  of  between  (20-30)%  based  on  post-test  measure¬ 
ments.  Therefore  we  could  not  use  the  simple  term  for  C  given  by 
Eqn. (3-5)  but  must  return  to  Eqn. (3-1)  for  our  analysis. 


*The  convention  of  Cq=  +1  for  a  black  object  is  widely  although  not 

universally  used;  some  use  C  =  -1. 

o 
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. 1  Selection  Rules 

The  Contrast  Attenuation  Technique  for  determining  dust  density  is 

J  m-  Koschmeider 1 s  Equation,  viz., 

C  =  C  exp(-BR) .  (3-6) 

o 

•i  the  MILL  RACE  experiment,  it  was  previously  shown  that  for  multiple 
i  -'M-wfs,  a  modified  form  of  Eqn .  (3-6)  produced  better  results,  e.g.. 


C. 

-  ~  -■  =  exp  [0(R.  -  R.)]  ,  i<j  .  (3-7) 

i 

We  will  examine  the  use  of  Eqn.  (3-7)  for  analyzing  the  MILL  RACE 
films  and  derive  a  few  simple  "Selection  Rules"  to  use  with  the  data  to 
minimize  processing  errors. 

Let  us  redefine  Eqn.  (3-7)  with  the  following  independent  variables: 


—  =  exp[3(R  -R  ) 3  (3-8) 

y  y  x 

where  x  and  y  are  two  contrasts  for  different  targets  on  the  same  frame  as 
measured  from  any  of  the  MILL  RACE  film. 

Solving  for  0,  we  find 


6  = 


R  -R 

y  x 


[in  x  -  in  y]. 


(3-9) 


Consider  a  small  incremental  change  in  3  and  its  relation  to  both  x  and  y , 


J  .e , 


A3  = 


R  -  R 

y  x 


fc-f] 


Divide  both  sides  by  Eqn.  (3-9)  to  obtain, 


l~Ax  -  A^l 

\3  =  Lx  yj  . 


43 

3 


(3-10) 


(3-11) 


in  x  -  in  y 


The  error  in  Ax  and  Ay  (AC^  and  AC.)  were  determined  by  the  film  processing 
and  the  resolution  of  the  densitometer  used  in  measuring  the  contrasts. 

Roth  of  these  factors  remain  constant  for  a  given  film  and  a  densitometry 
measurement  of  the  same  film  frame.  Therefore  we  can  equate  Ax  and  Ay;  thus, 


(“)  -  (-7±~) 

y 


AC, 


where  we  define  AC  =Ax  =  Ay=  Constant. 


(3-12) 
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The  term  on  the  left  of  Eqn.  (3-12)  is  the  incremental  error  in  the  deduced 
extinction  coefficient  due  to  the  measured  target  contrasts  for  two  tar¬ 
gets  x  and  y.  The  minimum  error  occurs  when 


a 


B 


M 

■ 


After  differentiating  the  right  hand  side  of  Eqn. (3-12) 
algebra)  we  find  the  following  expression: 


(3-13) 

(and  some  messy 


(3-14) 


Solving  Eqn.  (3-14),  we  find  that 


x 

y 

will  minimize  the  error  in 
Or 


3.01, 

A0 


(3-15) 


Ci  *  3.01  , 

C. 

J 

C.  =  (3.01)  C.  ,  (3-16) 

In  c.  =  in  C.  +  1.10 
i  J 

From  Eqn. (3-9),  we  have  for  minimum  error. 


=  (i-io) 

R.  -  R. 
.1  i 


(3-17) 


Although  Eqn.  (3-17)  gives  the  best  ratio  of  and  ,  it  does  not 

give  any  guidance  in  selecting  the  best  values  for  both  Ch  and  C.. .  We 
can  derive  this  by  looking  at  Eqn. (3-6)  , 

C  =  Co  exp(-BR) . 


or 


(3-18) 
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By  following  a  similar  argument  as  used  to  derive  Eqns.  (3-11)  to  (3-16), 
we  can  derive  the  following  equations: 


(3-19) 

(3-20) 


C  =  C  e  1  , 
o 

or 

C  =  ( .368) C  .  (3-21) 

o 

Returning  to  Eqn.  (3-12)  and  using  the  results  given  in  Eqn.(3-15) ,  we  can 

A3 

calculate  the  minimum  error  for  — g-  ■  -  ,  i  .e  . , 


AB  =  y  -  (3 .01)y  A  C 
3  (3.01)y  (1.10) 


or 

AB 

6 


-(0.607) 


A  C 
C 


For  the  MILL  RACE  data  analysis,  the  measurement  error  AC  was  approximately  0.02* 
From  Eqn. (3-21)  and  for  Cq-  0.7,  we  find  the  minimum  error  in  C  occurs  for 
C  =  .258  (3-22) 

Therefore,  we  find  the  minimum  error  in  the  deduced  attenuation  coefficient  to  be: 


or  about  5%. 

Based  on  the  results  shown  in  Eqns. (3-16)  and  (3-22),  we  have  derived  the 
following  guidelines  in  analyzing  the  MILL  RACE  data: 

1.  Apparent  Contrast  C 

0.1  <  C<  0.5 

2.  Apparent  Contrast  Ratio 

0.5<  6iC.  -  inC .  <  1.50. 

-  i  J~ 

During  data  analysis,  these  rules  were  observed  when  possible.  For  some 
cases,  the  rules  were  violated  because  of  a  lack  of  data;  these  cases  are 
identified. 


*  This  error  was  due  to  the  reading  resolution  of  the  densitometer  used  to 
measure  contrast  from  the  MILL  RACE  films. 
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DETERMINATION  OF  THE  CONTRAST  ATTENUATION  COEFFICIENT 
Equation  (3-7)  can  be  rewritten  as  follows: 


R.  -  R. 
J  i 


In 


C. 

i 

c. 

3 


(3-7) 


Applying  the  selection  rules  defined  in  the  ,'revious  section  to  the  con¬ 
trast  data,  we  deduced  values  for  8  at  both  the  10  psi  and  20  psi  station 
for  intervals  of  one  second  after  shock  arrival  until  40  seconds  and  at 
five  second  intervals  from  40  seconds  until  60  seconds.  The  results  are 
shown  in  Table  (3-2)  and  consist  of  an  average  and  standard  deviation  for 
B  when  more  than  one  pair  of  targets  meets  the  selection  rule  guidelines. 

All  of  the  data  shown  in  Table  (3-2)  were  obtained  from  the  Flight 
Research  35mm  data  films  by  using  the  diffuse  densitometry  measurements. 

In  addition,  we  made  measurements  on  several  selected  frames  of  the  GSAP 
16mm  data  film  at  the  10  psi  station  using  the  microdensitometer.  The 
results  are  shown  for  several  different  times  in  Table  (3-3).  Although  the 
agreement  is  not  exact,  the  differences  can  be  explained  by  considering 
the  methods  used  to  deduce  contrast  of  a  photo-target  with  the  diffuse 
densitometer  and  the  microdensitometer. 

Figure  (3-8)  depicts  a  photo-target  being  analyzed  by  both  instru¬ 
ments.  In  the  top  diagram,  we  see  that  the  densitometer  beam  covers  over 
one  half  the  image  of  the  photo-target  and  a  similar  area  for  the  back¬ 
ground  measurement.  Thus,  the  measured  film  contrast  is  deduced  by 
averaging  over  a  relatively  large  area.  Any  small  dust  inhomogeneities 
will  be  smoothed  out. 

On  the  other  hand,  the  bottom  half  of  Fig.  (3-8)  shows  the  size 
of  the  microdensitometer  beam  on  the  film.  Note  that  the  beam  size 
(.071mm  dia)  is  less  than  1/14  the  size  of  a  photo-target.  By  slightly 
displacing  the  position  of  the  beam,  it  was  possible  to  change  the 
measurement  of  the  film's  density  by  substantial  amounts.  In  particular, 
the  background  measurement  was  highly  sensitive  to  the  vertical  location 
of  the  densitometer's  beam  and  by  moving  the  beam  one  mm  on  the  16mm  frame, 
it  was  possible  to  change  the  contrast  reading  by  about  50%.  Obviously, 
with  that  sensitivity,  it  would  not  be  possible  to  expect  exact  agreement 
between  these  different  measurements. 


53 


As  an  additional  check  between  the  measurements  obtained  by  using 
either  the  diffuse  or  the  microdensitometer,  we  performed  the  following 
procedure.  Several  Flight  Research  film  frames  were  scanned  with  the  micro¬ 
densitometer  over  the  two  targets  and  backgrounds.  Since  we  previously 
obtained  data  from  these  films  with  the  diffuse  densitometer,  a  direct 
comparison  between  the  two  different  densitometers  on  the  same  film  was 
now  possible.  The  results  are  shown  in  Table  (3-4).  Again  we  see  that 
exact  agreement  does  not  occur  and  again  we  attribute  the  differences  to 
the  size  of  the  beam  used  to  scan  the  films  in  the  densitometer.  By  using 
the  difference  in  3  measured  by  the  microdensitometer  as  compared  to  the 
diffuse  densitometer,  we  calculated  the  required  change  in  background 
film  density  necessary  to  produce  this  difference;  the  results  ranged 
between  (5-8)%.  Since  a  typical  data  frame  had  a  background  film  density 
that  changed  by  over  a  factor  of  two  when  scanned  vertically  by  the  micro¬ 
densitometer,  a  change  of  even  10%  would  not  be  unreasonable  for  small 
changes  in  the  scanned  position. 


Table  3-2  (Continued) 


Time 

(sec) 


(standard 
S@10  psi  deviation) 
-1 


m 


35 

.045  (.004) 

.17  (.087) 

36 

.041  (.003) 

.16  (.114) 

37 

.047  (.005) 

— 

38 

.067  (.008) 

.086  — 

40 

.056  (.012) 

.11 

45 

.055  (.008) 

.080*  (.040) 

50 

.037  (.006) 

.024*  — 

55 

.048  — 

.014*  — 

60 

.030  — 

.015*  — 

(standard 
B@20  psi  deviation) 
-1 


m 


*  These  values  do  not  conform  to  the  selection  rules  given  in 
Section  3.3.1  and  therefore  used  with  discretion. 


Table  (3-3).  Attenuation  Coefficients  for  the  Flight  Research  35mm 
camera  and  the  GSAP  16mm  camera  at  10  psi  and  at  several  different 
times  after  shock  arrival. 


Table  (3-4).  Attenuation  coefficient  for  the  10  psi  Flight  Research 
camera  at  several  different  times  as  measured  by  the  diffuse  and  micro 
densitometer. 


Time 

8-  Diffuse  Densitometer 

3-  Microdensitometer 

(sec) 

(m  *) 

(m  *) 

11 

0.29 

0.43 

14 

0.30 

0.36 

28 

0.55 

0.47 

SECTION  4 

DETERMINATION  OF  DUST  DENSITIES 


In  Section  3  we  determined  values  for  the  attenuation  coefficient, 
3,  from  the  MILL  RACE  film  data;  these  data  are  shown  in  Table  3-1.  From 
Eqn.  (1-1)  (repeated  below)  in  order  to  determine  dust  densities,  we 

need  values  for  the  grain  density,  p  ,  and  the  term  enclosed  within  the 

S 

brackets . 


m 


4_ 

3 


°i6 


ra 

rna^. 

/  ~3 
a  . 
mm 


aJ  f (a) da 


max 

/  Q(a)a2f (a)da 
a  . 
min 


(l-D 


A  sample  of  dust  from  the  MILL  RACE  CAT  test  area  was  analyzed  by  members 

of  the  Atmospheric  Science  Laboratory  at  WSMR.  They  provided  a  value 

for  p  and  also  values  for  the  refractive  index  of  the  dust;  these  are 
8 

given  below. 


p  =  2.3  g/cm3  (4-1) 

8 

n  =  (1.54)-  i ( 2  x  10_3)  ,  (4-2) 

where  i  =  -y/-l  . 

In  Appendix  B,  the  rationale  for  choosing  Qext  instead  of  for 
use  in  Eqn.  (1-1)  is  given.  Fig.  (4-1)  shows  the  calculated  values  of 
Qext  that  were  based  on  Mie  theory,  the  refractive  index  given  in  Eqn.  (4-2), 
and  at  a  wavelength  of  A =0.55  ;  this  Q  was  used  for  calculating  the 
term  inside  the  brackets  of  Eqn. (1-1). 

The  expression  for  f(a)da  has  to  be  determined  from  the  MILL  RACE 
test  site  and  should  be  measured  at  each  interval  of  time  corresponding 
to  our  film  data.  A  separate  experiment,  fielded  by  SAI  with  the  assist¬ 
ance  of  the  U.S.  Army's  Atmospheric  Sciences  Lab,  was  tasked  to  make  these 
measurements.  A  separate  Report  in  these  Proceedings  describe  the  SAI 
experiment.  For  our  needs,  only  the  time-resolved  particle  size  distri¬ 
bution  function  was  necessary  to  deduce  dust  densities  using  Eqn.  (1-1) . 
However,  SAI  (Dr.  John  Cockayne)  volunteered  to  compute  the  entire 


c  n 


expression  within  the  brackets  using  our  values  for  Q(a)  and  a  previously 
supplied  computer  code  along  with  their  measured  values  of  f(a)da.  This 
was  done  and  the  results  are  shown  in  Table  (4-1)  for  the  10  psi  station. 

Given  the  data  of  Table  (4-1),  we  can  now  compute  absolute  dust 
density  vs  time  spectra  for  the  10  psi  station  using  Eqn.  (1-1).  The 
results  of  such  a  calculation  are  provided  in  Table  (4-2)  (where  the 
average  of  the  two  measurements  at  the  10  psi  station  was  used) .  The  cor¬ 
responding  values  of  dust  density  as  measured  by  SAI  are  shown  along  side 
our  CAT  deduced  data.  The  ratio  of  these  dust  densities  deduced  from  the 
two  different  methods  is  given  in  Table  (4-3).  Averaging  of  the  entries 
produces  the  following: 


0  10  psi 


p  (CAT) 
m _ 

p  (meas' 

_  m 


=  2.40  with  a  standard  deviation  of  0.61 


^ Average 


Referring  to  the  measured  data  of  Table  (4-1)  provided  by  SAI, 
the  two  measurement  instruments  at  the  10  psi  station  were  significantly 
different  in  their  readings;  over  the  same  time  interval  shown  in  Table  (4-2), 
the  two  measured  values  of  differed  by  an  average  factor  of  1.67  with  a 
standard  deviation  of  0.62,  i.e.. 


@  10  psi 


p  (meas  #1) 
m 

P  (meas  //2) 
m 


=  1.67  with  a  standard  deviation  of  0.62 


Average 


If  in  fact  the  dust  density  changed  by  a  factor  of  1.7  between  two 
measuring  instruments  separated  by  only  15  meters,  it  is  very  likely  that 
our  CAT-deduced  dust  density,  which  is  a  path-averaged  value  over  a  path- 
length  up  to  45  meters,  would  also  show  a  substantial  difference.  Thus, 
given  that  the  measured  values  for  the  dust  density  are  correct,  it  would 
not  be  unreasonable  to  expect  a  factor  of  2.4  between  the  CAT-deduced 
and  measured  dust  density. 
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Me  (4-2)  CAT  deduced  dust  densities  vs  time  and  measured  values  of  dust 
densities  vs  time  provided  by  SAI 


3 

Dust  Density  (10 _ g/cm  ) 


Time 

(sec) 


10  psi  station 

CAT  I  MEASURED 


5 

91 

7 

56 

7 

77 

7 

02 

6 

35 

4 

42 

3 

84 

4. 

82 

4 

35 

3. 

47 

3 

76 

3 

,32 

2. 

84 

2 

17 

2. 

29 

3 

31 

2. 

50 

2 

66 

1 

79 

1. 

47 

1. 

47 

0. 

82 

0 

55 

0. 

71 

0. 

87 

0. 

60 

0 

62 

0. 

56 

0. 

43 

0. 

37 

0. 

43 

0. 

66 

• 

0 

58 

0. 

53 

0. 

40 

0 

40 

0 

21 

Table  (4-3) 


Ratio  of  CAT-deduced  dust  density  C">  •’easured 
values  of  dust  density  vs  time 

p  (CAT)/u  (meas) 
m  ip 


Time 

(sec) 


I 


10  psj. 


static 


SECTION  5 
CONCLUSIONS 


The  results  of  the  MILL  RACE  validation  experiment  of  the  CAT 
method  for  deducing  dust  densities  indicate  that  this  method  produced 
value  greater  than  the  measured  dust  density  by  a  factor  of  2.4  ±  0.61 
By  using  reasonable,  conservative  arguments  we  believe  that  the  CAT 
method  produces  dust  densities  to  within  a  factor  of  about  two  of  the 
actual  path-averaged  values. 
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APPENDIX  A 

THEORETICAL  DERIVATION  OF  CONTRAST  ATTENUATION 
TECHNIQUE  EQUATIONS 

Consider  a  diffuse  target,  T,  located  at  distance  R  from  an  obser¬ 
vation  point,  0  (see  Fig.  A-l) .  Let  ^  be  the  solid  angle  at  0  subtended 
by  T  along  the  path  r. 

The  irradiance ,  E,  at  0  produced  by  the  target  is  given  by 


E(T)  .  m?  exp(-SexR> 


(A-l) 


where  is  the  extinction  coefficient  along  r  (assumed  constant)  and 
I  (T)  is  the  radiance  of  T.  For  visible  radiation  and  a  non-emitting 
target : 

I(T)  =  pH(R,t)  (A-2) 

where  p  is  the  reflectivity  of  T  and  H  (R,t)  is  the  visible  irradiance  of 
T  at  time  t  due  to  all  causes  (sun,  scattered  light,  etc.).  Thus, 

E(T)  =  pMliktLil  exP  (-eexR)  .  (A- 3 

A  second  source  of  visible  irradiation  of  0  arriving  within  Q  is  the 
so-called  "Airlight"  or  path  radiance  which  is  the  total  light  scattered 
into  the  field  of  view  of  0  along  r  from  all  sources  including  the  sun, 
sky  and  terrain.  The  airlight  irradiation  at  0  can  be  found  by  the 
following  expression: 

R 

E(A)  =  I  H(r,t)  8  exp  (-8  r)  dr,  (A-4 

q  sc  ex 

where  8  is  the  scattering  coefficient  towards  0.  For  the  case  when  8 

sc  ex 

H(r,t)  and  8gc  are  all  constant  along  r,  we  can  integrate  (A-4)  to  find: 
E(A)  -  H(t)  8sc  Cl-e*P<-SexW]. 


The  sum  of  (A-3)  and  (A-5)  describe  the  total  irradiance  at  0  arriving 
within  the  angle  Q: 
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(A-6) 


E(R)  =  E(T)  +  E(A)  =  pH(R,t)ft  exp(-8c.,  R) 

TT 

+  H(t)  8  [l-exp(-8  R)]. 

_ sc  ex 

8 

ex 

Note  that  H(R,t)  is  not  equal  to  H(t)  =  H(r,t  )  when  R  because  the 
target  itself  blocks  radiation  from  portions  of  the  surroundings  at  T. 

The  irradiance  at  0  produced  by  the  background  just  outside  of 
ft  is  found  from  Eqn.  (A-4)  with  the  limits  of  integration  extending  to 
infinity,  viz., 

CO 

E(B)  =  /  H(r,t)  8  exp  (-8  r)  dr 
o  sc  ex 


or  when  H(r,t,)  =  H(t) , 

E(B) 


H(T)B 

sc  . 


8 


ex 


(A-7) 


Equation  (A-7)  applies  when  the  background  radiance  is  not  due  to  any 
specific  feature  such  as  a  hill,  tree,  etc.,  but  due  only  to  airlight 
extending  into  the  sky  towards  the  horizon. 

Visibility  of  the  target  at  0  occurs  when  the  contrast  between  the 
target  and  its  background  exceeds  a  threshold  which  depends  on  the  detecting 
element.  Contrast,  C,  is  defined  by: 

r  _  E(R)  -  E(B)  ,  (A- 8) 

L  E(B) 


or, 


pH(R,t)ft 8  -  H(t)8 

—  ex  sc 

IT 


H(t)  sc 


exp(- 8  R)  . 
ex 


(A-9) 


The  inherent  contrast  of  a  target,  Cq,  is  found 

observer  at  0  approach  the  target,  i.e.,  letting  R-*-0. 

(A-l)  to  (A-8) ,  this  produces  the  following  expression: 

pH(R,t)ft  8  -  H(t)  6 

ex  sc 


by  letting  the 
Reviewing  Eqns. 

(A10) 


Co(T) 


H(t) 

SC 
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or  combining  Eqns.  (A-9)  and  (A-10) ,  we  find: 

C(t)  =  C  (t)  exp(-  3  R),  ( A- 1 1 ) 

o  ex 

which  is  commonly  referred  to  as "Koschmeider ' s  Law",  although  Duntley  and 

others  contributed  to  its  derivation. 

An  explicit  assumption  for  the  derivation  of  Eqn.  ( A— 11)  is  the 

uniformity  of  the  intervening  medium  between  the  target  and  the  observer. 

This  allows  the  use  of  a  single  extinction  and  scattering  coefficient, 

and  3gc,  to  be  used.  In  principle  one  could  derive  a  similar  expression 

with  non-uniform  3  and  3  :  in  this  case  one  would  have  to  know  their 
ex  sc 

properties  and  perform  several  additional  integrations  along  the  path  r. 

We  assume  that  an  average  value  of  3  and  3  exists  and  satisfies  Eqns. 

6X  SC 

(A-l)  to  ( A— 11)  . 

A  serious  limitation  on  using  Eqn.  ( A— 11)  to  find  3gxfor  a  time- 

varying  case  such  as  MILL  RACE  was  the  necessity  of  defining  both  C(t)  and 

C  (t)  at  the  same  time,  t.  For  cases  when  C  (t)  is  constant  or  well 
o  o 

defined,  one  can  use  pre-determined  measurements  and/or  calculations  for 

CQ(t)  along  with  the  measured  value  of  C(t)  and  thus  solve  Eqn.  (A-ll)  for 

3ex-  Unfortunately,  during  the  MILL  RACE  experiment  neither  case  applied, 

i.e.,  conditions  surrounding  the  target/camera  path  length  were  rapidly 

changing  in  time.  This  situation  affects  the  values  of  H(R,t) ,  H(t) ,  and 

also  values  for  3  and  3  .  An  alternative  method  can  be  used  if  the 

ex  sc 

following  conditions  existed: 

1.  Measurements  were  made  of  multiple  targets  at  the  same  time. 

2.  The  intervening  medium  was  uniform. 

3.  The  inherent  contrast  of  the  targets  were  all  equal. 

If  such  conditions  existed,  then  we  can  wri'te  an  Eqn.  (A-ll)  for 
each  independent  target,  i.e., 

C1 ( t)  =  CQ(t)exp(-3exR1) , 

c2(t)  =  CQ(t)exp(-3exR2) , 

etc.  Dividing  these  two  Eqns.,  we  can  write, 

C^t) 


r  ( 


=  exp(-BQjR,  "  RJ  ) 


( A— 12) 


s\-'' 


and  since  (R^  -  R2)  was  a  known  value  while  C^(t)  and  C2(t)  were  measured 

at  0,  we  can  solve  Eqn.  (A-12)  for  3  viz., 

ex 


I 


6ex(t)  - 


Cj(t ) 

c^TtjT 


R2  “  R1 


(A- 13) 


Of  course,  for  (n)  targets  we  can  solve  Eqn.  (A-12)  for  (n-1)  values  of 


3  viz . , 
ex. 


gex(t)  = 


[Ci<t>  1 

lCi+  (t)J  ,  i  =  1,2,...,  n-1. 


(A-14) 


R.  ,  ,  -  R. 
l+l  1 


Once  we  solve  for  3  ,  the  value  of  C  (t)  can  also  be  found  from  Eqn.  ( A— 11) 

ex  o 

By  iterating  this  technique  for  many  different  times  and  for  many  different 
targets,  it  becomes  possible  to  begin  a  systematic  study  of  the  time 
dependent  behavior  of  3 


■r.  v,v. 


APPENDIX  B 


MIE  SCATTERING  AND  EXTINCTION  CONSIDERATIONS 
Contrast  attenuation  of  visible  radiation  in  the  atmosphere 
occurs  according  to  Koschmeider' s  Eqn.: 

C  =  C  exp  (-  &R) ,  (B-l) 

K  O 

where  C  ,  C  are  the  apparent  contrast  at  distance  R  and  inherent  con- 
R  o 

trast,  respectively,  while  8  refers  to  the  attenuation  coefficient  for 

either  scattering  (8  or  extinction  (8  .  Since  the  basis  of  the 

scat  ext 

CAT  to  deduce  dust  densities  depends  upon  the  definition  of  8,  we  under¬ 
took  a  sensitivity  study  to  determine  the  best  choice  of  8  and  the  effect 
of  changing  its  value. 

A  brief  review  of  elementary  Mie  scattering  theory  is  necessary 

in  order  to  discuss  the  pertinent  factors.  Recall  that  Mie  theory  applies 

for  all  sizes  of  spherically  shaped ,  isotropic  particles  at  all  wavelengths. 

2TT  3. 

Mie  found  that  the  parameter  a=  — ^ -  (where  a  is  the  particle  radius  and 

A  the  wavelength)  determines  the  distribution  of  phase  over  the  particle; 
when  a<<l,  Mie  theory  reduces  to  Rayleigh  scattering  while  for  a>>10,  the 
principles  of  diffraction,  reflection  and  refraction  adequately  describe 
the  particle's  interaction  with  the  radiation.  In  the  midrange  of  rela¬ 
tive  particle  size,  however,  the  use  of  the  full  Mie  theory  becomes  manda¬ 
tory.  For  most  of  the  sizes  of  dust  particles  encountered  at  MILL  RACE, 
the  full  Mie  theory  applies. 

The  second  vital  parameter  used  for  Mie  scattering  calculations 
is  the  index  of  refraction  of  the  particle.  When  absorption  is  insigni¬ 
ficant,  only  scattering  occurs  and  the  refractive  index  is  expressed  as 
a  real  number,  n.  When  the  particle  is  absorbing  significantly,  the  re¬ 
fractive  index  must  be  expressed  by  a  complex  number,  i.e., 

m  =  n-in^  ,  (B-2) 

where  m  is  the  refractive  index  with  n  and  n^  the  real  and  imaginary  parts, 
respectively,  and  i  =n/TT  The  imaginary  part  is  related  to  the  absorption 
coefficient  of  the  particle. 

Referring  back  to  Koschmeider ' s  Eqn.,  we  can  relate  8gcat  to  the 
real  index  of  refraction  while  8  obtains  for  the  complex  case.  Once 


the  refractive  index  is  specified,  a  calculation  of  the  intensity  dis¬ 
tribution  function  for  radiation  interacting  with  the  particles  becomes 
tractable  through  the  use  of  the  Mie  theory.  The  results  are  usually 
expressed  in  terms  of  an  angular  cross  section,  a(0) ,  for  scattering  by 
the  particle,  and  defined  as  that  cross  section  of  an  incident  wave,  acted 
on  by  the  particle,  having  an  area  such  that  the  power  flowing  across  it 
is  equal  to  the  scattered  power  per  steradian  at  an  observation  angle  0. 

By  integrating  over  all  scattering  angles,  one  can  determine  the  total 
scattering  cross  section,  defined  as  that  cross  section  of  the  incident 
wave  acted  on  by  the  particle,  having  an  area  such  that  the  power  flowing 
across  it  is  equal  to  the  total  power  scattered  in  all  directions;  or 


4tt 

a  =  /  a  (0)du 

o 


(B-3) 


A  more  useful  parameter  than  a  is  the  efficiency  factor  Q,  which 
can  be  defined  for  both  scattering  and  extinction.  Qgc  is  defined  as  the 
ratio  of  the  scattering  cross  section  to  the  geometric  cross  section,  i.e. 


while  Q  _  is  defined  as 
ext 


Q  =  Q  +  Q  . , 
vext  sc  ab 


(B-4) 


(B-5) 


where  Q  ^  is  the  efficiency  factor  for  absorption.  Mie  theory  allows  one 
to  calculate  both  Q  and  Q  ;  their  difference  equals  Q  ,  .  We  do  not 

SC  6Xl  3D 

intend  to  discuss  further  the  derivation  of  these  parameters.  We  do 
intend  to  describe  the  differences  resulting  from  including  absorption 
in  the  CAT  theory. 


Figure  (B— 1)  shows  the  effect  of  including  absorption  in  the  cal¬ 
culation  of  Q  for  various  values  of  Q  ,  .  Note  that  as  absorption 
sc  ab 

increases  from  n^  <  10  4  to  n^  =  10  1 ,  the  limiting  value  of  Qgc  as  a 
becomes  very  large  is  between  1.00  and  1.11,  increases  to  1.2  when  n^  =  1, 
and  is  near  1.95  when  n^  =  10.  The  limiting  value  of  Q  approaches  two 
in  all  cases.  It  appears  that  the  asymptotic  values  of  both  Qgc  and  Qext 
approach  two  as  Q  becomes  larger,  while  for  very  small  values  of  Q  , 


the  limiting  value  for  Q  goes  to  one  while  Q  continues  towards  two. 

sc  ext 

Thus  for  very  large  particles  and/or  short  wavelength  cases,  we  need  to 
use  the  extinction  efficiency  factor  for  calculating  changes  in  contrast 
or  else  a  serious  error  (on  the  order  of  e2)  may  result. 

The  MILT.  RACE  data  pertains  to  dust  particles  with  radius  extend 
ing  upwards  of  30  microns.  For  visible  light,  this  produces  an  a^300 
which  satisfies  the  "very  large  a"  criteria.  Thus  we  need  to  compute  0 
for  this  case  to  determine  its  significance  in  calculating  the  contrast 
attenuation . 


Values  for  the  index  of  refraction  for  MILL  RACE  soil  were 
measured  by  WSMR  and  found  to  be : 

m  =  1.54  -  .0021 

A  calculation  was  made  for  both  Q  and  Q  at  260  different  values  oi 

sc  ext 

d  =  2a,  viz. , 


Ad  =  O.Oly  for  0.2<  d  <  1.0 

Ad  =  O.ly  for  1.0£  d  _<  10 

Ad  =  l.Ou  for  10  <  d  100, 


and  for  different  wavelengths, 


1  =  .40  .42  . . 68 , .  70 , 

A  =  .55y 

and  A  =  .6328  y 

(The  latter  wavelength  was  used  in  the  PMS  particle  sizing  instrument) . 

The  results  for  Q  and  Q  ^  at  A  =  .55  are  shown  in  Fig.  (B-2) .  As 
sc  ext 

expected,  the  asymptotic  values  diverge  towards  a  factor  of  two  difference, 
but  the  differences  at  small  particle  diameters  (d  _<  lOy)  remain  relatively 
close  together.  Fig. (B-3)  shows  the  normalized  difference  between  0gc  and  C 
The  importance  of  the  divergence  between  Qgc  and  Qext  depends  on 
the  relative  weighting  of  these  parameters  caused  by  the  particle  size 
distribution  function  used  to  calculate  8.  The  equation  for  8  is  the 
following : 


a 

max 

8  =  /  0  N  f(a)  da. 


(B-6) 


0  2  4  6  8  10  12  14  16  18  20  22  24  26  28 

a 


Fig.(B-l)  Efficiency  factor  for  scattering  as  a  function  of 
2ttb 

-  “  and  for  several  different  values  of  the  imaginary  com¬ 
ponent  of  the  refractive  index  n^.  The  real  component  of  the 
refractive  index  is  constant  and  equal  to  1.5.  (From  Appl. 
Opt.  5,  No.  2  February  1966  p.  281). 


where  N  is  the  total  particle  concentration,  f(a)da  is  the  fraction  of 

particles  having  radius  between  a  and  (a  +  da)  and  a  ,  a  .  are  the 

max  min 

maximum  and  minimum  size  particles  contributing  to  the  attenuation. 

We  can  rewrite  Eqn.  (B-6)  as: 


a 

max 

$  ~  ^  f  a2QNf(a)da.  (B-7) 

a  . 
mm 


We  see  that  the  effect  of  changing  Q  from  to  Qg  depends  on  the 

second  moment  of  the  distribution  function.  Likewise  the  determination 

of  p  ,  the  dust  mass  concentration  is  found  from  Eqn.  (B-8) : 
m  ^ 


a3f (a)da 


(B-8) 


P  =  T  P„ 

m  J  g 


max 

/  2 

a  QF(a)da 


where  0  is  the  grain  density,  and  the  other  terms  have  already  been 
S 

defined.  We  have  previously  shown  that  an  experimental  determination  of 


8  can  be  found,  e.g., 


(e)  . 


(B-9) 


R2  -  Ri 


so  that  the  derived  value  for  p  also  depends  on  the  value  of  Q  in  a  com- 

m 

Heated  way. 


APPENDIX  C 

SELECTED  DATA  ON  FILMS  USED  FOR 
THE  MILL  RACE  CAT  VALIDATION  TEST 


KODAK  HIGH  DEFINITION  AERIAL  FILM  3414 


(ESTAR  THIN  BASE) 


This  is  an  extremely  fine-grain,  slow  speed,  panchromatic, 
negative  aerial  camera  film. 

BASE:  A  2.5mil  (64  micrometer)  Estar  Thin  Base 

with  a  B2  (dyed-gel)  backing 


SENSITIVITY:  Panchromatic, 

AERIAL  EXPOSURE  INDEX: 
AERIAL  FILM  SPEED: 

RMS  GRANULARITY  VALUE: 
RESOLVING  POWER: 


with  extended  red  sensitivity 

Daylight  (no  filter)  -  2.5* 

Daylight  (no  filter)  -  8.0* 

9  (at  a  net  density  of  1.0)* 

Test-Object  Contrast  1000:1-630 

lines/mm* 

Test-Object  Contrast  1.6:1-250 

lines /mm* 


*These  data  are  based  on  development  in  Kodak  developer  D-19  at  66F  for 
8  minutes  in  a  sensitometric  processing  machine. 


KODAK  LINAGRAPH  SHELLBURST  FILM  2476 


(ESTAR  AH-BASE) 

This  is  a  medium  grain,  medium  speed,  panchromatic,  high 
resolution,  negative  film. 

BASE:  A  .004  inch  (102  micrometer)  ESTAR-AH  base  with  a 

fast-drying  backing 

SENSITIVITY:  Panchromatic  with  extended  red  sensitivity 
EXPOSURE  INDEX:  250* 

DIFFUSE  RMS  GRANULARITY  VALUE:  22* 

RESOLVING  POWER:  Test-Object  Contrast  1000:1-160 

lines /mm* 

Test-Object  Contrast  1.6:1-63 

lines /mm* 


*Based  on  development  in  Kodak  developer  D-19  at  68F  for  12  minutes  using 
continuous  agitation. 
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Naval  Surface  Weapons  Center 

ATTN:  Code  K06,  C.  Lyons 
ATTN:  Code  R15,  J.  Petes 
ATTN:  Code  F31 

Naval  Weapons  Evaluation  Facility 
ATTN:  Code  70,  R.  Tillery 
ATTN:  Code  70,  L.  Oliver 

Ofc  of  the  Deputy  Ch  of  Naval  Ops 

ATTN:  NOP  654,  Strat  Eval  &  Anal  Br 

Strat  Systems  Project  Office 
ATTN:  NSP-272 
ATTN:  NSP-273 
ATTN:  NSP-2722,  F.  Wimberly 
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DEPARTMENT  OF  THE  AIR  FORCE 

Aeronautical  Systems  Div 

ATTN:  ASD/ENFTV,  D.  Sorgen 
2  cy  ATTN:  ASD/ENFTV,  D.  Ward 

Air  Force  Geophysics  Lab 
ATTN:  LY,  C.  Touart 

Air  Force  Rocket  Propulsion  Lab 
ATTN:  LKLP,  G.  Beale 

Air  Force  Systems  Command 
ATTN:  XRTO 
ATTN  SDM 

Air  Force  Tech  Applications  Ctr 
ATTN:  TF 

Air  Force  Weapons  Lab 

ATTN:  NTYV,  A.  Sharp 
ATTN:  HO,  W.  Minge 
ATTN:  NTYV 
ATTN:  NTES 
ATTN:  SUL 
2  cy  ATTN:  NTO 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  FIMG 

ATTN:  FBAC,  D.  Resell us 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  MBC,  G.  Schmidt 

Air  University  Library 
ATTN:  AUL-LSE 

Arnold  Engrg  Dev  Ctr, 

ATTN:  AEDC,  DOFOV 

Air  Force  Systems  Command 

ATTN:  ENSN,  W.  Wilson 
ATTN:  SYDT 
ATTN:  EN 

ATTN:  HQ  Space  Dlv/RST 
ATTN:  ENMR 

ATTN:  HQ  Space  Div/RSS 
ATTN:  ENSN  Blankinship 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
ATTN  AFRDQI 
ATTN:  AFRD 

Foreign  Technology  Div 

ATTN:  SOBS,  J.  Pumphrey 
ATTN:  TQTD 
ATTN:  SDBG 

USAF  HDQ,  Directorate  of  Ops 
ATTN:  AFXOOTS 

Space  Division 

ATTN:  AFML 

Strat  Air  Command 


DEPARTMENT  OF  ENERGY 

Department  of  Energy 
ATTN:  OMA/RD&T 

OTHER  GOVERNMENT  AGENCY 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 


NATO  School,  SHAPE 

ATTN:  US  Doc  Officer 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  L-262,  J.  Knox 
ATTN:  L-125,  J.  Keller 
ATTN:  L-8,  R.  Andrews 

Los  Alamos  National  Lab 
ATTN:  R.  Thurston 
ATTN:  MS670,  T.  Scolman 
ATTN:  R.  Selden 
ATTN:  D.  Kerr 
ATTN:  MS  F668,  R.  Dingus 
ATTN:  J.  Hopkins 

Sandia  National  Lab 

ATTN:  Org  7112,  A.  Chabai 
ATTN:  M.  Cowan 

Sandia  National  Labs,  Livermore 
ATTN:  T.  Cook 

ATTN:  Library  &  Security  Class  Div 
ATTN:  H.  Norris 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

ACUREX  Corp 

ATTN:  C.  Nardo 
ATTN:  C.  Powars 

Aerojet  General  Corp 
ATTN:  R.  Steele 

Aerospace  Corp 

ATTN:  H.  Blase 
ATTN:  R.  Crolius 

Analytic  Services,  Inc,  ANSER 
ATTN:  J.  Selig 

APTEK,  Inc 

ATTN:  T.  Meagher 

AVCO  Systems  Division 
ATTN:  J.  Stevens 
ATTN:  W.  Reinecke 
ATTN:  W.  Brodlng 
ATTN:  J.  Gilmore 
ATTN:  P.  Grady 
ATTN:  A  .  Pall one 
ATTN :  Doc  Con 

Battelle  Memorial  Institute 
ATTN:  E.  Unger 
ATTN:  M.  Vanderlind 
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Boeing  Aerospace  Co 

ATTN:  M/S  13-13,  R.  Dyrdahl 

Boeing  Co 

ATTN:  M/S  41-52,  M.  Susman 

ATTN:  W.  Hanmon 

ATTN:  M/S  85/20,  E.  York 

ATTN:  R.  Holmes 

ATTN:  M/S  41/10,  J.  Avery 

Boeing  Military  Airplane  Co 

ATTN:  MS  75-74,  D.  Sawdy 

California  Research  &  Technology,  Inc 
ATTN:  K.  Kreyenhagen 
ATTN:  M.  Rosenblatt 

Cal  span  Corp 

ATTN:  M.  Holden 
ATTN:  M.  Dunn 

Dupont  Chemical  Corp 
ATTN:  F.  Bailey 

G.  B.  Laboratory,  Inc 

ATTN:  G.  Burghart 


Martin  Marietta  Denver  Aerospace 
ATTN:  E.  Strauss 

McDonnell  Douglas  Corp 
ATTN:  P.  Lewis  Jr 
ATTN:  L.  Cohen  MS  13-3 
ATTN-  R.  Reck 
ATTN:  G.  Johnson 
ATTN:  E.  Fitzgerald 
ATTN:  H.  Berkowitz 
ATTN:  D.  Dean 

McDonnell  Douglas  Corp 
ATTN:  M.  Potter 

National  Academy  of  Sciences 

ATTN:  National  Materials  Advisory  Board 
ATTN:  D.  Groves 

Northrop  Corp 

ATTN:  B.  Butler 
ATTN:  J.  Eucs 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 
ATTN:  G.  Lang 


General  Electric  Co 
ATTN:  P.  Cline 
ATTN:  B.  Maguire 

General  Research  Corp 
ATTN:  J.  Mate 

H-Tech  Labs,  Inc 

ATTN:  B.  Hartenbaum 


Harold  Rosenbaum  Associates,  Inc 
ATTN:  G.  Weber 


Hercules,  Inc 

ATTN:  P.  McAllister 


Institute  for  Defense  Analyses 
ATTN:  Classified  Library 

Kaman  AviDyne 

ATTN:  R.  Ruetenik 
ATTN:  S.  Criscione 
ATTN:  N.  Hobbs 


Pan  Am  World  Service  Inc 

ATTN:  AEDC/Library  Doc,  TRF 

PDA  Engineering 

ATTN:  J.  Schutzler 
ATTN:  M.  Sherman 
ATTN:  J.  Dunn 
ATTN:  J.  McDonald 

Physics  International  Co 
ATTN:  J.  Shea 

R&D  Associates 

ATTN:  W.  Graham 
ATTN:  P.  Rausch 
ATTN:  J.  Carpenter 
ATTN:  F.  Field 

Rand  Corp 

ATTN:  R.  Rapp 

Rockwell  International  C  -p 
ATTN:  G.  Perroue 


Kaman  Sciences  Corp 

Rockwell  International  Corp 

ATTN: 

J.  Harper 

ATTN: 

R.  Hemann 

ATTN: 

J.  Hoffman 

ATTN: 

F.  Shelton 

S-CUBED 

ATTN: 

J.  Keith 

ATTN: 

R.  Duff 

ATTN: 

G.  Gurtman 

Kaman  Science  Corp 

ATTN: 

D.  Sachs 

Science  Applications,  Inc 

ATTN: 

W.  Yengst 

Kaman  Tempo 

ATTN: 

W.  Plows 

ATTN: 

B.  Gambill 

ATTN: 

C.  Lee 

ATTN: 

DAS I AC 

ATTN: 

J.  Warner 

ATTN: 

J.  Manship 

Lockheed  Missiles  S  Space  Co,  Inc 

ATTN: 

J.  Stoddard 

ATTN: 

F.  Borgardt 

Science  Applications,  Inc 

Lockheed  Mlsslfes  &  Space  Co,  Inc 

ATTN: 

W.  Lay  son 

ATTN: 

R.  Walz 

ATTN: 

J.  Cockayne 

-FiD-A134  723 


UNCLASSIFIED 


VALIDATION  OF  THE  CONTRAST  ATTENUATION  TECHNIQUE  (CAT) 
FOR  DEDUCING  DUST  <U>  TELEDVNE  BROWN  ENGINEERING 
HUNTSVILLE  AL  M  L  PRICE  28  FEB  82  SD81-DNA2574 
DNA-TR-81-81  DNA001-81-C-8129  F/G  14/5 
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Science  Applications,  Inc 

TRW  Electronics 

&  Defense  Sector 

ATTN:  A.  Martellucci 

ATTN 

A. 

Zimnerman 

ATTN 

T. 

Williams 

Southern  Research  Institute 

ATTN 

N. 

Lipner 

ATTN:  C.  Pears 

ATTN 

P. 

Brandt 

ATTN 

R. 

Bacharach 

SRI  International 

ATTN 

M. 

King 

ATTN:  D.  Curran 

ATTN 

W. 

Wood 

ATTN:  G.  Abrahamson 

ATTN 

M. 

Seizew 

ATTN:  H.  lindberg 

ATTN 

R. 

Plebuch 

ATTN 

A. 

Ambrosio 

System  Planning  Corp 

ATTN 

D. 

Baer 

ATTN:  J.  Jones 

ATTN 

T. 

Mazzola 

Terra  Tek,  Inc 

TRW  Electronics 

&  Defense  Sector 

ATTN:  S.  Green 

ATTN 

D. 

Kennedy 

ATTN 

V. 

Blankinship 

Teledyne  Brown  Engrg 

ATTN 

D. 

Glenn 

4  cy  ATTN:  M.  Price 

ATTN 

W. 

Pol ich 

ATTN 

E. 

Wong 

Thiokol  Corp 

ATTN 

P. 

Dai 

ATTN:  J.  Hinchman 

ATTN 

L. 

Berger 

ATTN:  W.  Shoun 

ATTN 

N. 

Guiles 

Toyon  Research  Corp 

ATTN:  J.  Cunningham 

ATTN:  B.  Gragg 

ATTN 

E. 

Allen 
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